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SUMMARY 


The  successful  design  of  a  propeller-driven  V/STOL 
aircraft  requires  the  precise  calculation  and  associated 
optimization  of  the  propeller  performance  in  static  opera¬ 
tion.  A  general  theory  for  performance  prediction  has  been 
formulated  based  on  a  continuous  vortex  representation  along 
the  lines  of  a  classical  lifting-line  model. 

As  opposed  to  forward  flight,  the  deformation  of  the 
trailing  vortex  wake  is  appreciable  just  behind  the  propel¬ 
ler,  and  its  determination  constitutes  the  heart  of  the 
static  thrust  problem.  Emphasis  has  been  placed  on  deter¬ 
mining  a  satisfactory  force-free  approximation  to  the  ef¬ 
fective  "pitch”  of  the  trailing  vortex  sheets  with  the 
contraction  pattern  fixed  according  to  a  heavily  loaded 
actuator  disk  theory. 

Numerical  techniques  and  associated  computer  programs 
have  been  developed  to  calculate  not  only  the  inflow  to  the 
propeller  but  also  the  velocity  induced  along  arbitrarily 
described  deformed  trailing  vortex  sheets.  In  the  numeri¬ 
cal  calculations  for  a  series  of  specific  propellers  and 
representative  blade  loadings,  iterations  were  made  on  the 
effective  "pitch"  variations  as  well  as  on  the  inflow. 
Comparisons  of  predicted  and  measured  performance  were  not 
coupletely  satisfactory  so  far  as  the  "pitch"  iterations 
were  pursued,  but  indications  were  that  the  amount  of  con¬ 
traction  should  be  iterated  as  well. 

Techniques  for  gathering  and  reducing  instantaneous 
hot-wire  measurements  of  the  wake  velocity  were  developed. 
Measurements  were  carried  out  for  two  propellers.  Non- 
repeatable  data  were  obtained  in  large  regions  behind  the 
propeller  tip  in  both  cases.  Theoretical  calculations 
could  not  be  iterated  successfully  in  these  cases,  so  sig¬ 
nificant  theoretical  and  experimental  comparisons  could  not 
be  made. 

Detailed  observations  have  been  made  on  the  character 
of  the  numerical  computations,  and  certain  generalizations 
have  been  made  which  lead  to  computational  simplifications. 
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INTRODUCTION 


High  performance  in  the  static,  or  hovor,  operating 
condition  is  a  key  design  condition  for  V/STOL  propellers 
because  the  total  thrust  must  exceed  the  gross  weight  of 
the  aircraft  for  VTOL  capability.  Furthermore,  the  payload 
of  these  vehicles  is  typically  201  of  the  gross  weight,  so 
that  any  improvement  in  the  static  performance  results  in 
a  large  percentage  of  increase  in  payload  capacity. 

The  design  of  high-hove’ ing- performance  propellers 
has  been  hindered  by  the  lack  of  reliable  theoretical 
methods  for  predicting  the  static  performance.  The  prin¬ 
cipal  difficulty  in  developing  a  static  theory  is  in 
finding  the  force-free  location  of  the  trailing  vortex 
system.  This  is  in  contrast  to  the  forward  flight  theory, 
whore  the  location  of  the  trailing  vortex  system  is 
determined  primarily  by  the  free- stream  ’'elocity  and 
propeller  rotation.  In  forward  flight,  the  axial,  radial, 
and  tangential  deformations  of  the  sheets  due  to  the  induced 
velocity  do  not  occur  sufficiently  close  to  the  propeller 
nor  in  great  enough  magnitude  to  have  a  major  effect  on  the 
performance.  In  static  operation,  however,  these  deforma¬ 
tions  occur  immediately  behind  the  propeller  and  are  of  such 
magnitude  lhat  they  must  be  accurately  represented.  In 
fact,  the  axial  induced  velocity,  which  leads  only  to  a 
perturbation  of  the  pitch  of  'he  regular  helical  trailing 
vortex  sheets  in  forward  flight,  determines  the  entire 
variation  of  the  effective  "pitch"  in  static  operation. 

In  initial  studies  of  the  static  thrust  problem,  we 
used  a  simple  wake  model  in  which  the  essential  deformation 
was  discontinuous;  see  Refs.  25,  7,  8  and  9.  The  results 
were  generally  unsatisfactory,  bu*  they  did  indicate  that 
the  smooth,  continuous  deformation  of  the  trailing  vortex 
system  must  be  adequately  represented.  We  began  inves¬ 
tigation  of  a  continuous  deformation  model  under  U.S.  Army 
AVLABS  Contract  DA  44- 177-AMC- 165(T)  and  reported  the 
results  in  Refs.  10  and  11.  Although  a  so-called  refined 
wake  hypothesis  developed  in  that  study  as  an  approximation 
to  the  force-free  condition  showed  promise  of  reasonable 
agreement  between  performance  predictions  and  test  data,  it 
was  concluded  that  the  hypothesized  wake  velocities  should 
be  checked  by  actual  computations  as  well  as  by  test 
measurements . 


In  the  present  investigation,  we  have  proceeded  in 
these  directions.  In  conjunction  with  Canadair,  Ltd.,  and 
Dr.  I.  S.  Gartshore  and  Mr.  D,  C.  Gilmore  of  McGill 
University,  the  instantaneous  velocity  components  in  the 
wakes  of  two  propellers  have  been  measured  (Ref.  14).  More¬ 
over,  numerical  techniques  have  been  developed  to  compute 
the  induced  velocity  components  along  the  trailing  vortex 
sheets.  These  techniques  have  been  used  to  continue  the 
development  of  a  force-free  approximation  that  is  satisfac¬ 
tory  for  reliable  performance  predictions. 

The  theoretical  formulation  presented  in  Chapter  1  is 
nearly  identical  to  that  of  our  earlier  investigation  (Refs. 
10  and  11),  since  the  research  described  in  this  report  is 
a  logical  development  of  the  continuous  deformation  model 
of  tne  earlier  work.  The  details  are  repeated  here  to  insure 
that  this  report  is  self-contained  with  all  the  assumptions 
and  approximations  clearly  evident. 


CHAPTER  1 


THEORETICAL  FORMULATION 


Basic  Equations 

We  consider  a  propeller  rotating  at  a  constant  angular 
speed  n  in  a  uniform,  inviscid  fluid  which  has  a  density  p 
and  is  at  rest  at  infinity.  The  blade  thickness-to-chord 
and  chord- to- radius  ratios  are  assumed  to  be  sufficiently 
small  that  the  classical  lifting-line  formulation  should  be 
a  reasonable  approximation.  We  can  then  represent  the  pro¬ 
peller  blades  by  radial  bound  vortex  lines  accompanied  by 
a  force-free  system  of  trailing  vortex  sheets. 

A  propeller- f ixed  cylindrical  coordinate  system  (x,  r, 
0)  is  chosen  with  the  axis  of  rotation  as  the  x-axis;  see 
Figure  1.  The  N  blades  are  located  in  the  plane  x  ■  0  and 
are  arranged  symmetrically  such  that 

9p  ■  2w(p-l)/N  p  -  1,  2,  ...  ,  N  (1) 

i.e.,  the  p  ■  1  blade  coincides  wit.i  0-0. 

The  induced  velocity  anywhere  in  the  propeller  field 
is  determined  by  integration  of  the  Biot-Savart  law  over 
all  the  vortex  elements  comprising  the  system  in  terms  of 
their  strength,  orientation,  and  distance  to  the  field 
point  of  interest,  and  may  be  expressed  in  the  form 


u 


v 
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p-i  Rh 


[  rCrp)  UB 

[r(V  VB 
£r(rp)  W„ 


dr(rp) 
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dr(rp) 
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where 


u  is  thu  axial  component  of  the  velocity  induced  at  a 
field  point  (x,  r,  0)  and  taken  positive  in  the 
positive  x-direction, 

v  is  the  corresponding  radial  component,  positive  in 
the  positive  r-direction,  and 

w  is  the  corresponding  tangential  component,  positive 
in  the  positive  e-direction; 

is  the  radius  of  the  propeller  hub; 

Rp  is  the  radius  of  the  propeller  tip; 

r  is  the  bound  blade  circulation  strength  at  (0,  rD, 

ep) ;  p 

dr  is  the  strength  of  the  vortex  sheet  element  trailing 
from  CO,  rp,  ep) ; 

rp  is  the  radial  variable  of  integration  along  the  p**1 
p  propeller  blade; 

UB  is  the  influence  function  for  the  axial  component  of 
the  velocity  induced  at  (x,  r,  0)  by  a  bound  radial 
vortex  element  of  unit  strength  and  unit  length  at 

(°»  *p,  ep); 

Vg  is  the  corresponding  function  for  the  radial  compo¬ 
nent  ,  and 

Wg  is  the  corresponding  function  for  the  tangential 
component ;  and 

U.p  is  the  influence  function  for  the  axis1  component 
of  the  velocity  induced  at  (x,  r,  0)  by  an  arbi¬ 
trarily  deformed  vortex  sheet  element  of  unit 
strength  and  semi- infinite  length  trailing  from 
(0 ,  rp ,  0p) , 

Vj  is  the  corresponding  function  for  the  radial 
component,  and 

WT  is  the  corresponding  function  for  the  tangential 
component . 
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The  influence  functions  can  be  derived  from  first 
principles  (Ref.  27,  pp.  1.10  -  1.12,  for  example),  or  they 
can  be  obtained  from  previous  results  (Ref.  23,  pp.  13  -  16). 
They  are  summarized  in  Table  I,  where 


x 


V 


is  the  axial  coordinate  of  any  point  on  the 
trajectory  of  the  vortex  sheet  element  trailing 
from  (0,  rp,  9p) , 

is  the  corresponding  radial  coordinate,  and 
is  the  corresponding  azimuthal  coordinate; 
is  the  value  of  u  at  (xv,  rv,  ev) , 
is  the  corresponding  value  of  v,  and 
is  the  corresponding  value  of  w;  and 
is  a  time  parameter  to  be  defined. 


When  a  point  on  a  lifting  line  is  chosen  as  a  field 
point,  the  contribution  from  the  bound  vortex  for  this 
blade  is  irrelevant  and  so  discarded,  as  in  wing  theory 
(Ref.  15,  pp.  131  -  133). 


The  contribution  from  the  bound  vortices  to  the  induced 
velocity  is  simple  to  calculate.  On  the  other  hand,  the 
contribution  to  the  induced  velocity  from  the  elements  of 
the  trailing  vortex  sheets  is  very  complicated.  Not  only 
do  the  influence  functions  require  an  integration  over  the 
length  of  the  element,  but,  what  is  even  worse,  the  location 
or  "trajectory"  of  the  element  is  not  fixed.  Rather,  these 
elements  must  drift  force  free  by  aligning  themselves  with 
the  streamlines  of  the  flow,  which  in  differential  form 
(Ref.  27,  pp.  1.1  -  1.2,  for  example)  are 


dx  .  dr  ,  de 

u  v  fi+  (w/r)  (5) 


These  equations  can  be  integrated  in  terms  of  the  parametric 
time  t  that  it  takes  a  fluid  particle  to  move  along  this 
trajectory,  giving 


x 


V 


(6) 
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SUMMARY  OF  INFLUENCE  FUNCTIONS 


(7) 


t 

0V  -  0p  ♦  JU  ♦/  (wv/ry)  dx 


The  integrands  in  Eqs.  (6)  through  (8)  are  functions  of  the 
coordinates  along  the  trajectory,  and  the  coordinates  them¬ 
selves  are  functions  of  t. 


For  convenience,  we  now  introduce  a  set 
such  that  Eqs.  (2)  through  (4)  become  simply 

of  operators 

u  -  0u(r, -dr/dr;  u,v,w) 

(9) 

v  -  Ov(r, -dr/dr;  u,v,w) 

(10) 

w  ■  0w(r, -dr/dr;  u,v,w) 

(ID 

where 

0U  is  the  Biot-Savart  operator  for  the  axial 
induced  velocity  component, 

Ov  is  the  corresponding  operator  for  the  radial 
component ,  and 

Ow  is  the  corresponding  operator  for  the  tangential 
component. 

Note  carefully  that  the  operators  operate  on  the  bound  and 
trailing  vortex  strengths  with  the  velocity  field  itself 
appearing  parametrically.  These  equations  constitute  a  set 
of  simultaneous,  nonlinear,  singular  integral  equations  over 
the  domain  of  the  blades  and  the  trailing  vortex  sheets. 

This  set  is  complete  if  r  is  specified. 

If  r  is  not  specified,  another  equation  is  required  to 
relate  the  circulation  strength  and  velocity  field  through 
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the  propeller  geometry.  Consider  the  force  and  velocity 
diagram  for  a  typical  section  of  the  blade  at  the  radial 
station  (0,  Tp,  0p),  Figure  2,  where 

Wp  <s  the  total  local  velocity  seen  by  the  blade 
v  section, 

Up  is  the  local  axial  inflow  [i.e.,  the  axial  compo- 
F  nent  of  the  induced  velocity  at  (0,  Tp,  8p)], 

v  is  the  corresponding  radial  inflow  (not  shown,  but 
F  perpendicular  to  the  plane  of  the  section)  ,  and 

Wp  is  the  corresponding  tangential  inflow; 

a  is  the  local  blade  angle  of  attack; 

B  is  the  local  blade  pitch  setting; 

+p  is  the  angle  between  nrp  and  Wp; 

y  is  the  tangent  angle  of  the  sectional  lift- 

drag  polar; 

dD  is  the  elemental  profile  drag  on  the  blade 
section;  and 

dL  is  the  elemental  lift  on  the  blade  section. 

From  the  definition  of  the  sectional  lift  coefficient 
in  terms  of  the  local  blade  chord  b,  namely 

CL  2(dL/drp)/pWp2  b  (12) 

the  Kutta- Joukowski  formula  (Ref.  27,  pp.  5.21  -  5.22,  for 
example)  can  be  written  as 


r  -  b  Wp  Cl/2  (13) 

where  the  lift  coefficient  is  assumed  to  be  known  (Ref.  6, 
for  example)  as  a  function  of  a  for  a  given  section  at  ap¬ 
proximately  the  local  Mach  and  Reynolds  numbers.  This  is 
the  equation  needed  in  conjunction  with  Eqs.  (9)  through 
(11)  to  complete  the  set  when  the  propeller  geometry  is 
specified. 
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Figure  2.  Velocity  and  Force  Diagram  at  Typical 


In  either  case,  the  circulation  or  the  propeller 
geometry  specified,  the  prediction  of  the  overall  perfor¬ 
mance  can  be  computed  once  the  inflow  has  been  calculated. 
In  particular,  if  we  resolve  the  resultant  of  dL  and  dD  as 
indicated  by  the  dotted  lines  in  Figure  2  and  integrate 
along  the  blade  for  the  N  blades,  the  total  thrust  T 
and  power  P  in  nondimensional  form  are 
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where  the  thrust  and  power  coefficients  of 

the  propeller 

are  defined  by 

CT 

«2T/4Pn2Rp4  - 

T/pn2Dp4 

(16) 

CP 

n3P/4pfl3Rp5  ■ 

P/pn3Dp5 

(17) 

respectively,  with 

n  »  2n  n  and 

Rp  -  Dp/2. 

From  the  thrust  and  power  coefficients  we  can  calculate 
the  figure  of  merit  F/M  (Ref.  21,  pp.  353  -  355),  given  in 
percent,  by 


F/M  =  lOO^f  2CT3/irCp2  ■  79.8  C^,2/C?  (18) 

To  conclude  this  section,  it  is  of  interest  to  discuss 
the  nondimensionalization  of  the  basic  equations  of  the 
static  thrust  problem.  As  we  have  implied  by  the  form  of 
Eqs.  (16)  and  (17),  the  characteristic  time  is  l/il  and  the 
characteristic  length  is  Rp.  If  the  equations  are  nondi- 
mensionalized  with  respect  to  these  quantities,  complete 
similarity  would  exist  were  it  not  for  the  sectional  lift 
and  drag  data,  which  in  general  must  include  the  effects 
of  Mach  number  and  Reynolds  number,  as  noted  earlier. 
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Outline  of  Solution 


The  set  of  equations  described  above  cannot  be  solved 
analytically,  so  approximate  numerical  techniques  must  be 
employed.  Iteration  is  perhaps  the  most  straightforward 
approach  and  the  one  which  we  have  chosen. 


First  we  rewrite  Eqs.  (9) 
but  equivalent,  form;  namely, 


through  (11)  in  a  different, 


u  ■  u 

v  ■  v 

w  ■  w 


Ku[u 

*  ou(r, 

...)] 

(19) 

Kv[v 

-  ov(r, 

...)] 

(20) 

Kw[w 

-  ow(r, 

...)] 

(21) 

where 

K  is  a  factor  chosen  as  necessary  to  achieve 
convergence  of  the  axial  component  of  the 
induced  velocity  field, 

Ky  is  the  corresponding  factor  for  the  radial 
component ,  and 

Ky  is  the  corresponding  factor  for  the  tangential 
component. 

These  iteration  factors  may  be  constant  or  may  depend  on  any 
of  the  variables  of  the  problem. 

With  the  propeller  geometry  specified,  iteration  begins 
by  finding  a  zeroth,  or  initial,  approximation  to  the  inflow 
and  thus  the  blade  circulation  through  Eq .  (13),  the  Kutta- 
Joukowski  formula.  The  first  approximation  to  the  inflow 
and  wake  velocity  field  follows  by  operating  with  the  Biot- 
Savart  operator  on  the  zeroth  approximation  to  the  blade 
circulation  and  some  guess  at  the  induced  velocity  field  in 
the  propeller  wake.  These  results,  along  with  the  iteration 
factors  and  the  zeroth  approximations  to  the  inflow  and  the 
wake  velocity  field,  are  then  substituted  into  the  right- 
hand  sides  of  Eqs.  (19)  through  (21)  to  give  the  first  ap¬ 
proximations  to  the  inflow  and  the  wake  velocity  field, 
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respectively.  The  iteration  continues  by  operating,  as 
before,  on  the  first  approximations  to  determine  a  second, 
and  so  forth.  Once  convergence  within  a  prescribed  accuracy 
is  reached,  the  propeller  performance  can  be  calculated. 

On  the  other  hand,  for  the  case  with  the  blade  circu¬ 
lation  specified,  iteration  would  proceed  similarly  except 
that  the  circulation  would  not  change  at  each  step.  Upon 
convergence  in  this  case,  the  blades  could  be  designed 
consistent  with  the  Kutta - Joukowski  formula,  and  the  perfor¬ 
mance  could  then  be  calculated. 

The  iteration  approach  necessitates,  at  each  step  of 
the  iteration,  calculation  of  the  induced  velocity  field  at 
a  sufficient  number  of  points  to  define  accurately  the  local 
streamlines  in  the  propeller  wake.  Unfortunately,  the 
number  of  calculations  required  for  adequate  representation 
of  the  complete  distributed  vortex  sheets  is  very  large. 

As  a  result,  we  looked  for  s:me  way  to  reduce  the  number  of 
calculations  but  still  satisfy  the  force-free  condition  of 
the  wake,  at  least  in  an  approximate  sense.  This  search  for 
a  simplified,  yet  accurate,  iterative  approach  to  the  pre¬ 
diction  of  static  propeller  performance  has  been  the  goal 
of  both  the  earlier  and  the  present  investigations. 

The  question  of  convergence  and  the  proper  choice  of 
the  iteration  factors  in  Eqs.  (19)  through  (21)  have  proved 
to  be  important  and  difficult.  They  will  be  treated  in 
depth  in  our  discussion  jf  the  numerical  results  in  Chapters 
4  through  6. 


Continuous  Deformation  Model 

The  trailing  vortex  sheets  of  a  propeller,  in  general, 
stretch  axially,  together  with  a  radial  contraction  and  a 
tangential  distortion.  The  edges  of  the  sheets  are  also 
locally  unstable  and  tend  to  roll  up  as  for  a  finite  wing. 

While  these  deformations  are  not  of  primary  signifi¬ 
cance  for  the  practical  determination  of  aircraft  propel¬ 
ler  performance  for  the  usual  forward  flight  operating 
conditions  (Ref.  2),  they  are  the  heart  of  the  static  prob¬ 
lem.  Without  a  free  stream,  not  only  do  the  deformations 
occur  much  closer  to  the  propeller  and  so  exert  more  influ¬ 
ence,  but  they  are  also  much  larger.  In  fact,  the  axial 
elongation  which  is  only  a  perturbation  of  the  basic  helical 
sheets  in  forward  flight  now  determines  the  total  variation 
of  the  effective  "pitch"  of  the  sheets. 
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At  the  beginning  of  our  earlier  investigation  (Refs. 

10  and  11)  we  decided,  on  the  basis  of  previous  studies 
(Refs.  7  through  9  and  25),  that  there  was  a  clear  lack  of 
fundamental  knowledge  on  the  nature  of  the  roll-up  process 
for  propellers.  Moreover,  we  had  found  that  it  was  important 
to  represent  especially  well  the  axial  and  radial  deformation 
of  the  distributed  trailing  vortex  sheets  immediately  behind 
the  propeller.  Therefore,  we  decided  at  that  time  to  con¬ 
centrate  on  the  continuous  axial  stretching  and  radial  defor¬ 
mation  of  the  trailing  vortex  sheets  or,  in  other  words,  a 
"Continuous  Deformation  Model".  In  the  present  investigation 
we  have  continued  to  develop  this  model  with  principal  em¬ 
phasis  on  improvement  of  the  force- free  approximation,  espe¬ 
cially  with  regard  to  the  axial  deformation.  This  develop¬ 
ment  proceeded  throughout  the  numerical  computations  of  this 
study  and  is  described  in  Chapters  4  through  6,  including 
comparisons  with  test  data  wherever  possible. 


Implementation 

In  order  to  carry  out  even  approximate  iterative 
solutions  to  the  *et  of  equations  that  we  have  formulated, 
it  is  necessary  to  use  digital  computation  techniques, 
especially  for  the  Biot-Savart  operators  in  Eqs.  (2)  through 
(4).  Accordingly,  we  have  developed  a  very  general  computer 
program  (CODEFVEL) ;  see  Chapter  2  of  Ref.  10.  This  program 
evaluates  the  Biot-Savart  operators  for  field  points  (x,  r, 

0)  located  on  the  lifting  line  or  anywhere  else  except  on 
the  trailing  vortex  sheets  proper.  CODEFVEL  also  contains 
many  options  for  treating  wake  models  other  than  the  Con¬ 
tinuous  Deformation  Model.  All  of  this  generality  leads, 
unfortunately,  to  extremely  time-consuming  and  complicated 
data  preparation  (see  Appendix  I  of  Ref.  10). 

In  the  present  investigation  we  have  developed  two 
separate,  simplified  programs  that  are  specialized  to 
compute  either  the  inflow  at  field  points  along  the  lifting 
line  only  (INFLOWPT),  or  the  velocity  induced  at  field 
points  downstream  along  the  trailing  vortex  sheets  only 
(SHEETPT) .  Both  of  these  programs  are  based  upon  C0DEFVEL1 
but  they  do  contain  many  improvements  that  were  indicated  • 
by  our  experience.  INFLOWPT,  SHEETPT,  and  the  associated 
numerical  analysis  for  each  are  described  in  full  in  Chapters 
2  and  3,  respectively. 

Once  the  inflow  has  been  computed,  the  propeller  circu¬ 
lation  distribution  and  performance  can  be  found  from  Eqs. 
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CHAPTER  2 


METHOD  OP  INFLOW  AND  PERFORMANCE  COMPUTATION 
Description  of  Computer  Programs 

Program  INFLOWPT  was  developed  from  CODEFVEL,  as 
mentioned  in  Chapter  1,  to  provide  a  simpler,  quicker  means 
of  evaluating  the  inflow  Biot-Savart  operators  for  arbitrary 
axial,  radial,  and  tangential  deformation  of  the  trailing 
vortex  sheets  in  the  Continuous  Deformation  Model.  Just  as 
we  had  developed  a  short  program  (VORTTRAJ)  from  CODEFVEL 
for  computing  only  selected  trajectories  of  the  trailing 
vortex  elements,  we  found  it  convenient  to  abstract  an 
analogous  short  program  (MRCHINFO)  from  INFLOWPT.  MRCHINFO 
has  the  additional  option,  if  desired,  of  computing  the 
inflow  influence  functions  of  the  trailing  elements  as  well 
as  their  trajectories. 

Program  PROPERFM  was  written  as  an  independent  program 
for  evaluating  the  propeller  blade  circulation  and  perfor¬ 
mance.  Although  provisions  were  made  for  its  eventual 
inclusion  in  INFLOWPT  as  a  series  of  subroutines,  this  was 
not  practical  for  the  numerical  computations  that  were 
actually  performed. 

The  programs  were  all  written  in  FORTRAN  63  for  the 
CDC  1604  computer.  The  particular  machine  that  we  used  was 
located  at  the  Cornell  University  Computing  Center.  A 
duplicate  FORTRAN  63  deck  of  each  program  has  been  cut  and 
supplied  to  USAAVLABS  along  with  a  set  of  complete  listings. 
Brief  operating  instructions  are  provided  in  Appendix  I. 


Representation  of  Circulation  Distribution 

In  general,  the  blade  circulation  distribution  r  is 
tabulated  numerically  at  several  radial  stations,  but  for 
the  subsequent  integrations  in  Eqs.  (2)  through  (4)  we  need 
r  in  analytic  fora.  A  suitable  form  is  a  Glauert-type 
series  (Ref.  15,  pp.  138  -  139),  or 


r  •  nR^  g£  sint  *p  (22) 
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rp  =  1/2 (Rp+Rh)  -  l/2(Rp-Rh)  cos 


P 


(23) 


where 

Gfc  is  the  A**1  nondiraensional  Glauert  coefficient,  and 

"**P  is  the  Glauert  variable  which  runs  from  0  to  it 
with  its  center  midway  between  the  hub  and  tip. 

This  form  possesses  a  square-root  behavior  at  both  the 
tip  and  the  hub.  For  the  tip,  this  is  the  proper  behavior 
required  for  a  lifting- line  formulation  (Ref.  5  pp.  171  -  174, 
for  ex^mole) .  For  the  hub,  it  is  approximately  correct  when 
the  blade  necks  down  rapidly  at  the  shank,  as  on  some  propel¬ 
lers.  On  the  Other  hand,  when  the  hub  is  like  the  root  of 
a  wing,  the  slope  of  the  circulation  distribution  should  be 
zero  from  a  rigorous  viewpoint.  This  could  have  been  achieved 
exactly  by  an  additional  condition  on  the  determination  of 
the  Glauert  coefficients,  but  practically  the  contribution 
of  the  hub  to  the  overall  performance  is  so  small  that  zero 
slope  can  be  achieved  with  or  without  this  condition. 

The  Glauert  coefficients  are  evaluated  by  the  method 
of  least  squares,  which  minimizes  the  error  at  the  specified 
data  points  (Ref.  26,  pp.  363  -  370).  INFLOWPT  provides  for 
the  numerical  values  of  r  to  be  given  at  up  to  100  points, 
with  L  as  high  as  25.  Usually  L  -  10  was  the  best  compro¬ 
mise  between  the  accuracy  of  the  r  fit  and  the  introduction 
of  extraneous  fluctuations  into  the  shape  of  its  derivative. 
This  point  will  be  discussed  at  greater  length  in  the  section 
on  checks  of  the  program  and  also  in  the  section  on  the 
calculations  for  the  12.6°  pitch  setting  of  the  65  AF  propel¬ 
ler  in  Chapter  6. 


Influence  Functions  for  a  Trailing  Vortex  Element 

As  pointed  out  earlier,  we  must  in  general  calculate 
the  trajectories  of  the  elements  of  the  trailing  vortex 
sheets  in  order  to  evaluate  the  contribution  by  these  sheets 
to  the  induced  velocity.  The  method  described  here  is  used 
in  INFLOWPT  for  each  element  of  a  distributed,  continuous 
representation  of  the  trailing  vortex  sheets.  However,  the 
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method,  and  the  influence  functions  computed  by  it,  would 
be  equally  applicable  to  a  discrete  vortex  representation 
of  the  trailing  sheets. 


Physically,  most  of  the  wake  deformation  occurs  just 
downstream  of  the  propeller,  say,  within  a  distance  of  a 
blade  radius  or  so.  Therefore,  numerical  integration  of 
Bqs.  (6)  through  (8)  is  necessary  only  in  this  region,  since 
farther  downstream  the  trajectories  may  be  approximated  as 
regular  helices  with  suitable  pitch. 

To  carry  out  the  trajectory  integrations,  a  prescribed 
velocity  field  along  the  trailing  vortex  sheets,  either  from 
a  guess  or  from  the  preceding  step  of  the  iteration,  has  to 
be  stored  in  the  computer  as  input  data.  Though  the  real 
flow  field  is  not  ax i symmetric,  we  need  not  consider  the  e 
variation  because  we  are  always  moving  along  trajectories 
that  are  the  same  for  all  blades  at  a  given  radius.  The 
stored  velocity  field  is  specified  at  a  maximum  of  30  values 
each  of  x  and  r  and  3-point  Lagrangian  interpolation,  first 
in  r  and  then  in  x,  carried  out  for  the  other  values  in 
between.  We  found  that  for  each  of  12  radial  stations,  9 
axial  stations  were  adequate  to  represent  the  velocity  field 
between  the  propeller  plane  and  1.5  propeller  radii  downstream. 


With  the  velocity  field  stored,  the  trajectory  is  evalu¬ 
ated  in  a  step-by-step  marching  scheme  using  a  two-term 
constant  slope  method.  The  first  step  off  the  blade  is 
made  by  multiplying  the  inflow  components  by  a  small,  but 
arbitrary,  time  increment  At  to  find  the  new  location  of  the 
vortex  sheet  element.  The  second  and  successive  steps  are 
taken  by  interpolating  to  determine  the  velocity  at  this  new 
location,  multiplying  it  by  2  At,  and  adding  the  result  to 
the  coordinate  of  the  previous  location.  With  this  scheme, 
the  error  is  of  the  order  of  the  cube  of  At  as  opposed  to 
the  square  of  At,  which  an  ordinary  two- term  Taylor  series 
would  afford.  The  increment  At  is  constant  as  the  march 
proceeds,  but  the  increment  may  be  changed  to  new  values  at 
up  to  nine  points  downstream  to  speed  the  march  by  increasing 
and  decreasing  the  At  increment  where  allowable.  At  each 
change  in  increment  size,  the  two- term  constant  slope  march 
begins  again  as  it  does  from  the  blade. 

Simultaneously  with  this  march,  it  is  convenient  to 
perform  the  integrations  over  t  for  evaluation  of  the  con¬ 
tribution  of  this  deformed  trailing  vortex  element  to  the 
influence  functions;  see  Table  I.  The  integrands  are  cal¬ 
culated  at  each  new  location  in  the  march  for  the  element 
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from  each  blade,  are  summed  over  all  the  blades,  and  are 
integrated  by  the  trapezoidal  rule.  Running  totals  of  the 
influence  functions  are  carried  as  the  march  proceeds. 


The  trajectories  and  influence  functions  from  this 
scheme  were  compared,  by  means  of  MRCHINFO,  with  corre¬ 
sponding  results  from  CODEFVEL.  After  a  series  of  inves¬ 
tigations,  we  settled,  for  a  four-bladed  propeller,  on  45 
marching  steps  at  OAt  ■  0.01,  63  at  0.05,  and  the  remainder 
at  0.10.  For  a  three-bladed  propeller,  we  took  55  steps  at 
0.01,  40  at  0.05,  and  the  remainder  at  0.10.  The  march  along 
each  trajectory  automatically  terminates  when  either  all  of 
the  specified  steps  have  been  taken  or  the  downstream  cutoff 
limit  of  the  input  velocity  field  has  been  reached,  which¬ 
ever  occurs  first.  In  practice,  we  always  specified  suffi¬ 
cient  steps  at  fiAt  -  0.10,  say,  500,  so  that  each  march 
terminated  at  a  specified  cutoff  point  1.5  propeller  radii 
downstream.  This  means  of  ending  the  march  provides  a  great 
time  savings  over  CODEFVEL,  in  which  all  elements  took  the 
same  number  of  steps.  This  number  was  dependent  on  insuring 
that  the  "slowest"  element  reached  the  cutoff. 

When  the  termination  point  is  reached,  each  trajectory 
is  issumed  to  be  of  regular  helical  shape  with  pitch  deter¬ 
mined  from  the  axial  and  tangential  velocity  components  at 
that  point.  The  computation  of  the  influence  functions  can 
then  be  completed  by  approximating  the  contributions  of  this 
regular  helical  part  that  trails  downstream  from  the  termi¬ 
nation  point  to  infinity.  This  approximation  consists,  as 
in  CODEFVEL, of  computing  the  e-average  of  the  contributions 
using  results  of  the  theory  of  the  Generalized  Actuator  Disk 
(Ref.  17,  pp.  13  -  15).  The  calculation  involves  simply  the 
evaluation  of  some  well-defined  function;  namely,  the 
Legendre  functions  of  the  second  kind  and  degree  plus  and 
minus  one-half  and  the  Heuman  lambda  function.  We  chose  to 
evaluate  the  Legendre  functions  by  relating  them  to  the 
complete  elliptic  integrals  of  the  first  and  second  kinds 
(Ref.  3,  p.  249),  which  are  computed  using  a  rational  ap¬ 
proximation  (Ref.  1,  pp.  591  -  592).  For  the  Heuman  lambda 
function,  we  took  an  exact  series  expansion  (Ref.  3,  p.  300) 
and  related  it,  also,  to  the  complete  elliptic  integrals  of 
the  first  and  second  kind.  These  integrals  are  calculated 
as  before.  Where  applicable,  the  addition  formula  for  the 
Heuman  lambda  function  (Ref.  3,  p.  36)  is  employed  to  speed 
convergence  of  the  series. 
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Inflow  Induced  by  the  Trailing  Vortex  Sheets 


Once  the  influence  functions  have  been  computed, 
integration  of  the  second  terms  in  Eqs.  (2)  through  (4)  must 
be  carried  out.  We  have  chosen  to  perform  these  integrations 
for  a  continuous,  distributed  trailing  vortex  sheet  repre¬ 
sentation.  Unfortunately,  the  singularities  of  the  integrand 
complicate  the  picture,  and  special  treatment  is  required  in 
order  to  carry  out  these  integrations  quickly  and  accurately. 

We  first  see  that  the  square-root  behavior  of  Eq.  (22) 
makes  the  strength  of  the  elements  of  the  trailing  vortex 
sheets  infinite  at  the  hub  and  the  tip.  This  is  no  problem, 
though,  if  we  carry  out  the  integrations  in  terms  of  the 
Glauert  variable  instead  of  the  original  radial  variable. 

A  second  and  more  serious  difficulty  is  the  existence 
of  singularities  which  arise  in  the  influence  functions 
because  each  inflow  point  coincides  with  an  element  of  the 
trailing  vortex  system.  The  nature  of  these  singularities 
can  be  determined  by  direct  investigation  of  the  integrands 
of  the  influence  functions  in  Table  I,  but  the  actual 
strength  of  the  singularities  is  dependent  upon  the  particu¬ 
lar  path  along  which  the  inflow  point  is  approached.  This 
path  must  be  determined  by  the  location  of  the  elements  of 
the  trailing  vortex  sheet  immediately  adjacent  to  the  inflow 
point,  since  it  is  the  induction  of  these  elements  that  gives 
rise  to  the  singularity.  Thus,  for  inflow  points,  this  path 
is  the  radial  coordinate  along  the  lifting  line.  This  radial 
direction  is  tangent  to  the  trailing  vortex  eheet  at  the 
lifting  line  but,  in  general,  is  not  normal  to  the  individual 
trailing  elements  due  to  their  contraction.  However,  the 
resulting  singularity  will  be  nonintegrable  over  the  radius 
unless  the  radial  inflow  is  zero  there,  so  that  the  trailing 
elements  are  normal  to  the  lifting  line.  This  is  the  same 
nonintegrable  singularity  that  has  arisen  in  attempts  to 
extend  lifting- line  theory  to  swept-back  wings  (Refs.  13,  28, 
and  20,  for  example).  We  handled  this  difficulty  by  allowing 
the  first  few  steps  in  the  marching  scheme  to  be  taken  with 
zero  radial  velocity  and  constant  axial  and  tangential  veloc¬ 
ity,  i.e.,  as  regular  helical  elements,  after  which  the 
trailing  elements  deform  in  the  desired  fashion.  Generally 
we  have  taken  two  such  steps.  We  will  discuss  this  point 
further  in  the  section  on  the  calculations  for  the  12.6° 
pitch  setting  of  the  65  AF  propeller  in  Chapter  6. 

With  elimination  of  this  nonintegrability ,  we  are  left 
with  two  types  of  singularities,  Cauchy  and  logarithmic,  as 
shown  elsewhere  (Ref.  24,  pp.  14  -  16).  To  handle  the 
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integration  over  these  singularities,  we  break  the  results 
into  "singular"  and  "regular"  parts.  Once  the  singularity 
strengths  have  been  determined  analytically  from  the  basic 
equations  by  approaching  the  inflow  point  along  the  proper 
radial  path,  we  can  define 


flr+w  1 


4* [u^+ (flr+w) 1/2  rp-r 


(nr+w)3 _ 

8irr[u2+(nr+w)2]3/2 


in 


(24) 


(25) 


u 

nr+w 


(26) 


where  u  and  w  are  the  input  inflow  at  r.  These  singular 
parts  then  can  be  integrated  analytically  in  the  usual 
principal  value  sense  with  the  derivative  of  r  found  by 
differentiation  of  Eq.  (22). 

The  regular  parts  for  each  element  of  the  trailing 
vortex  sheets  are  determined  numerically  by  subtracting  the 
singular  parts  from  the  summation  over  all  the  blades  of  the 
influence  functions.  The  radial  location  of  the  elements 
is  chosen  automatically  in  INFLOWPT  by  specifying  as  input 
an  integer  number  of  basic  divisions  of  the  0  to  *  interval 
of  integration  over  the  Glauert  variable.  Since  the  major 
part  of  the  induction  is  by  the  outboard  trailing  elements, 
the  accuracy  of  integration  is  improved  by  automatically 
halving  the  basic  division  size  between  2w/3  and  * ;  i.e., 
outboard  of  approximately  the  0.8  radial  station.  Vortex 
elements  are  also  automatically  added  close  to  either  side 
of  each  inflow  point  to  provide  better  definition  of  the 
regular  parts.  These  elements  are  used  for  that  inflow 
point  only,  and  any  basic  elements  falling  nearer  the  inflow 
point  are  excluded  in  the  interest  of  numerical  accuracy. 

The  regular  parts  of  the  total  integrand  can  be  found  at  any 
desired  station  by  multiplying  the  regular  parts  of  the 
influence  functions  there,  as  found  from  the  computed  values 
by  3-point  Lagrangian  actual  integration  over  the  Glauert 
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variable  in  carried  out  by  a  generalized  Gaussian  integration 
scheme.  In  this  scheme ,  a  10-point  Gaussian  quadrature  is 
first  applied  to  the  total  interval  of  integration  The 
total  interval  is  in  turn  subdivided  into  two  intervals,  the 
10-point  Gaussian  quadrature  is  applied  to  each  of  these 
intervals,  and  the  results  are  added  to  get  the  answer. 

Next,  the  total  interval  is  subdivided  into  three  intervals, 
and  so  on.  At  each  step  after  the  second,  the  answer  is 
compared  with  the  two  values  which  immediately  precede  it, 
and  the  subdivision  is  continued  until  these  three  values 
are  within  a  desired  accuracy. 

Program  INFLOWPT  can  accommodate  up  to  200  trailing 
vortex  elements  and  can  compute  the  inflow  at  up  to  40 
points.  Generally  we  considered  8  inflow  points  with  basic 
elements  every  tt/6,  giving  a  total  of  25  trailing  elements. 
The  accuracy  specified  for  the  Gaussian  integration  was 
-  0.00005  of  the  tip  rotational  speed. 


Calculation  of  Propeller  Performance 

Program  PROPERFM  is  used  to  compute  the  propeller  blade 
circulation  distribution  and  performance  by  means  of  the 
equations  in  Figure  2  and  Eqs.  (13)  through  (15)  and  (18). 
The  essential  inputs  are  the  blade  geometrical  parameters, 
the  two-dimensional  sectional  lift  and  drag  data  at  ap¬ 
proximately  the  local  Mach  and  Reynolds  numbers,  and  the 
inflow.  These  data  are  read  in  at  each  desired  blade 
section  with  the  lift  coefficients  over  a  range  of  angles 
of  attack,  as  well  as  the  drag  coefficients  either  over  a 
range  of  lift  coefficients  or  over  a  range  of  angles  of 
attack  given  in  tabular  form.  It  is  necessary  to  prepare 
these  sectional  data  in  advance  by  hand.  We  used  the  basic 
Curtiss-Wright  data  (Ref.  6)  extended  to  high  angles  of 
attack  by  the  methods  of  Ref.  12.  Other  data  could  be  used 
equally  well,  but  small  changes  would  have  to  be  made  in  the 
program.  The  integrations  in  Eqs.  (14)  aid  (15)  are  carried 
out  over  the  radial  variable  by  means  of  the  generalized 
Gaussian  integration  scheme  described  above. 


Checks  of  Computer  Programs 

Many  checks  were  made  in  order  to  gain  confidence  in 
programs  INFLOWPT,  MRCHINFO,  and  PROPERFM,  as  well  as  to 
determine  the  best  choice  of  the  input  parameters  for 
insuring  quick,  yet  accurate,  results. 
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Initially  we  checked  MRCHINFO  and  1NFLOWPT  by  dupli¬ 
cating  cases  that  we  had  run  by  CODEFVEL  for  both  three- 
and  four-bladed  propellers.  Successful  completion  of  these 
checks  gave  us  confidence  in  the  compatibility  of  the  new 
programs  with  the  earlier  one. 

An  independent  check  was  afforded  by  comparing  our 
results  with  existing  ones  for  a  particular  marine  propel¬ 
ler  in  cruise.  We  had  made  calculations  earlier  of  the 
induced  velocity  at  field  points  external  to  the  blades 
and  regular  helical  wake  of  this  propeller  for  the  Naval 
Ship  Research  and  Development  Center  using  CODEFVEL.  Excel¬ 
lent  agreement  was  found  by  S.  B.  Denny  (Ref.  4)  between 
our  results  and  comparable  ones  computed  with  a  program 
developed  by  J.  E.  Kerwin  (Ref.  18).  For  completeness  we 
also  checked  both  CODEFVEL  and  INFLOWPT  by  computing  the 
inflow  to  this  propeller.  The  agreement  of  the  axial  inflow 
with  that  of  Kerwin  was  everywhere  within  2. SI,  with  the 
worst  discrepancies  near  the  hub  and  tip.  These  were  traced 
to  differences  in  the  fit  of  r  and  were  an  early  indication 
of  this  sensitivity,  which  we  will  discuss  below. 

Besides  this  comparison  of  inflow  calculations,  we 
also  found  excellent  agreement  between  our  computed 
influence  functions  for  regular  helical  elements  and  tabu¬ 
lated  results  by  W.  B.  Morgan  (Ref.  22).  We  concluded  that 
our  computational  methods  are  indeed  adequate  in  the  case 
of  propellers  in  cruise  with  regular  helical  trailing  vortex 
sheets . 

Once  the  computational  scheme  and  computer  programs 
had  been  checked  in  this  fashion  against  our  previous 
results  as  well  as  against  independent  ones,  we  examined  the 
effect  on  the  inflow  of  varying  certain  input  parameters. 

By  means  of  numerical  experimentation  with  MRCHINFO  and 
INFLOWPT,  we  determined  the  choices  that  we  have  described 
in  the  previous  section?*  for  the  numbers  of  marching  steps 
and  their  associated  increment  sizes,  the  number  of  elements 
needed  to  represent  the  distributed  trailing  vortex  sheets, 
the  number  of  points  at  which  the  input  velocity  field  must 
be  specified,  and  the  number  of  Glauert  coefficients  in  the 
blade  circulation  fit. 

The  most  important  of  these  is  the  number  of  Glauert 
coefficients  for  fitting  r.  Numerical  experiments  with  the 
marine  propeller,  for  which  we  did  not  have  r  data  outboard 
of  the  0.9  radial  station,  indicated  that  small  percentage 
changes  in  r  between  the  0.9  station  and  the  tip  can  lead 
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to  much  larger  percentage  changes  In  the  singular  part  of 
the  axial  inflow  in  the  same  region.  This  occurs  because 
the  singular  part  of  the  axial  inflow  near  the  tip,  and  hub 
as  well,  is  very  sensitive  to  the  higher  order  Glauert  co¬ 
efficients,  being  roughly  proportional  near  the  hub  and  tip 
to  the  summation  to  L  of  the  products  of  the  coefficients 
times  the  square  of  their  order  i.  In  contrast,  r  is  propor¬ 
tional  to  the  summation  of  the  coefficients  alone,  and  its 
derivative  is  proportional  to  the  summation  of  the  coeffi¬ 
cients  times  their  order.  Thus,  small  changes  in  the  r  fit 
near  the  hub  or  tip  can  change  the  higher  order  coefficients 
enough  that  the  changes  in  r  become  amplified  in  the  singu¬ 
lar  part  of  the  inflow.  Unfortunately,  this  occurs  near 
the  tip  where  the  accuracy  of  the  inflow  is  most  important. 

We  checked  PROPERFM  by  comparing  the  predicted  perfor¬ 
mance  of  the  Curtiss-Wright  3 (13168A10P3)  propeller  with  the 
hand  calculations  report  earlier  (Ref.  10,  p.  52).  The  dif¬ 
ferences  were  on  the  order  of  1  to  21  in  thrust  and  power 
coefficients  and  about  1/21  in  figure  of  merit.  They  arise 
due  to  the  change  from  the  trapezoidal  integration  by  hand 
to  the  generalized  Gaussian  scheme  in  PROPERFM. 
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METHOD  OF  VORTEX  SHEET  VELOCITY  COMPUTATION 


Description  of  Computer  Program 

Program  SHEETPT  was  developed  from  both  CODEFVEL  and 
INFLOWPT  to  provide  a  means  of  computing  the  Biot-Savart 
operators  for  field  points  located  on  trailing  vortex  sheets 
which  may  have  arbitrary  axial,  radial,  and  tangential  defor¬ 
mation.  The  basic  method  of  computation  is  identical  to 
INFLOWPT  in  most  respects.  Exceptions  are  that  the  actual 
location  of  the  field  points  on  the  sheets  must  be  determined, 
the  singular  parts  appropriate  to  these  points  evaluated,  and 
the  velocity  induced  by  the  bound  blade  vortices  calculated. 
All  of  these  differences  are  described  in  this  chapter. 

This  program,  too,  was  written  in  FORTRAN  63  for  the 
CDC  Computer,  and  a  duplicate  deck  and  listing  have  been 
supplied  to  USAAVLABS.  Brief  operating  instructions  are 
provided  in  Appendix  I. 


Determination  of  Sheet  Point  Locations 


The  first  task  in  computing  the  velocity  along  the 
trailing  vortex  sheets  is  to  determine  the  actual  location 
(x,  r,  0)  of  the  field  points.  It  is  desirable  for  our 
purposes  to  calculate  the  velocity  at  several  axial  positions 
x  along  the  trajectories  from  each  of  several  different 
radii,  rQ,  say,  along  the  lifting  line.  Therefore,  to  deter¬ 
mine  r  and  0  for  each  x  and  r«,  we  must  march  along  each  of 
these  elements  exactly  as  we  ao  in  the  calculation  of  the 
influence  functions.  When  the  step  is  taken  that  goes  beyond 
the  first  specified  x,  the  march  stops,  and  a  3-point 
Lagrangian  interpolation  with  this  point  and  the  two  preceding 
ones  is  used  to  determine  r  and  0  at  this  x.  Next,  the 
velocity  components  are  found  at  this  x  as  at  any  point  on 
the  sheets  by  the  basic  interpolation  in  the  input  velocity 
field.  Finally,  a  numerical  differentiation  technique  based 
on  3-point  Lagrangian  interpolation  is  used  to  find  the  time 
rate  of  change  along  the  trajectory  of  the  velocity  compo¬ 
nents  at  the  field  point.  These  will  be  used  later  in 
computing  the  logarithmically  singular  parts  of  the  influence 
functions.  The  march  then  resumes  until  r  and  0  have  been 
determined  similarly  for  all  specified  x  locations. 


Program  SHEBTPT  provides  for  following  elements  down¬ 
stream  to  as  many  as  10  axial  stations  from  as  many  as  20 
radial  stations.  Generally,  we  considered  four  typical 
radial  stations  and  four  axial  stations.  It  turned  out  to 
be  more  efficient  in  computing  costs  to  make  several  runs 
for  each  case  with  each  run  treating  one  or  two  axial 
locations . 


Influence  Functions  for  a  Trailing  Vortex  Element 

The  marching  scheme  for  determining  the  trajectories 
and  influence  functions  in  SHEETPT  is  nearly  identical  to 
the  one  in  INFLOWPT.  There  is  one  significant  difference, 
though,  to  insure  accurate  integration  of  the  influence 
functions  when  the  trajectories  are  close  to  the  field 
points.  The  march  proceeds  as  usual  until  it  reaches  a 
specified  axial  distance  from  each  field  point  x  location. 
There  the  march  automatically  changes  to  a  smaller  At 
increment,  which  is  used  until  the  march  has  reached  an  equal 
distance  beyond  the  field  point.  Then  the  march  resumes  and 
continues  as  usual.  This  occurs  around  each  x  location. 
Outside  of  these  intervals,  the  step  numbers  and  increment 
sizes  are  specified  as  in  INFLOWPT.  The  march  terminates 
identically  to  INFLOWPT,  and  the  regular  helical  contribu¬ 
tions  to  the  influence  functions  are  computed  in  the  same 
way . 


After  some  numerical  experimentation,  we  settled  on  a 
series  of  step  and  increment  sizes  which  we  used  for  both 
three-  and  four-bladed  propellers.  In  runs  with  axial 
locations  x  of  0.05  and  0.2  propeller  radius,  we  used  five 
marching  steps  at  flat  ■  0.01  and  500  at  0.05  initially  off 
the  lifting  line,  500  at  0.05  between  the  0.05  and  0.2 
stations,  and  70  at  0.0S  and  500  at  0.10  downstream  of  the 
0.2  station.  Within  a  Ax  of  0.043  propeller  radius  on  either 
side  of  the  field  points,  OAt  -  0.01  was  used  and  the  march 
was  terminated  at  1.7  propeller  radii  downstream;  i.e.,  1.5 
radii  beyond  the  0.2  station  to  be  consistent  with  INFLOWPT 
practice.  In  runs  where  only  the  axial  location  x  of  0.5 
propeller  radius  was  considered,  we  took  five  steps  at  0.01 
and  500  at  0.05  initially  and  70  at  0.05  and  500  at  0.10 
downstream  of  0.5  radius.  Within  0.060  propeller  radius 
on  each  side  of  0.5,  QAt  ■  0.01  was  used  and  the  march  was 
terminated  at  2.0  propeller  radii  downstream.  Finally,  for 
axial  locations  x  of  1.0  propeller  radius,  we  used  five  steps 
at  0.01,  20  at  0.10,  and  500  at  0.05  initially  with  70  at 
0.05  and  500  at  0.10  downstream  of  the  1.0  radius.  Within 
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0.060  on  either  side  of  1.0  propeller  radius,  we  used  QAt  - 
0.01  and  terminated  the  march  at  2.5  propeller  radii  down¬ 
stream. 

| 

Vortex  Sheet  Self- Induced  Velocity 

i 

Integration  of  the  second  terms  of  Eqs.  (2)  through  (4) 
is  carried  out  for  the  sheet  points  in  exactly  the  same 
fashion  as  for  the  inflow.  Singularities  again  occur  when 
the  sheet  point  and  an  element  of  the  trailing  vortex  sheets 
coincide.  Here,  too,  we  can  separate  the  integrand  into 
singular  and  regular  parts  so  that  the  integration  can  be 
carried  out  exactly  as  for  the  inflow. 

The  nature  of  the  singularities  can  be  found  by  direct 
investigation  of  the  integrands  in  Table  I,  and  their 
strengths  of  the  singularities  can  be  found  be  approaching 
the  sheet  point  (x,  r,  0)  along  the  path  appropriate  to  the 
adjacent  elements  of  the  trailing  vortex  sheet.  To  find 
this  path,  we  examine  the  curvilinear  line  which  is  formed 
by  the  intersection  of  the  trailing  vortex  sheet  with  the 
plane  passing  through  the  sheet  point  normal  to  the  x-axis; 
see  Figure  3.  The  appropriate  path  is  then  along  this 
curvilinear  line  or,  in  the  limit  as  the  point  (x,  r,  0)  is 
approached  and  the  singularity  strengths  are  determined, 
along  the  tangent  to  the  line  at  the  point. 

!  1 

This  path  is  not  normal  to  the  trajectory  through  the 
sheet  point  because  the  radial  velocity  is  generally  not 
zero.  However,  this  is  not  a  problem  in  this  case,  because 
the  trajectory  extends  upstream  of  the  point  as  well  as 
downstream,  and  so  the  singularity  is  integrable.  It  is 
precisely  this  feature  that  makes  J.  Weissinger's  (Ref.  31) 
approximation  of  placing  the  lifting  line  at  the  one-quarter- 
chord  line  and  satisfying  the  downwash  at  the  three-quarter- 
chord  line  possible  for  finite  swept  wings.  The  singularities 
are  integrable  at  the  three-quarter-chord  line,  whereas  they 
are  not  at  the  one-quarter-chord  line  where  the  trailing 
vortex  sheet  originates. 

Numerical  determination  of  this  path  of  approach 
requires  computation  of  the  angle  y  ,  which  the  curvilinear 
intersection  line  makes  with  the  local  radius  through  the 
sheet  point.  This  is  done  by  carrying  out  marches  to 
determine  the  r  and  0  of  two  closely  adjacent  points  on  the 
intersection  line.  Once  these  are  found,  y  is  computed  by 
numerical  differentiation  of  a  3-point  Lagrangian  inter¬ 
polation  through  these  two  added  points  and  the  sheet  point. 
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Intersection  of  Trailing  Vortex  Sheet  With 


A  suitable  choice  for  the  location  of  these  added  trajec¬ 
tories  is  coincident  with  those  particular  vortex  elements 
which  are  added  to  provide  better  definition  of  the  regular 
parts,  just  as  in  INFLOWPT. 

The  strengths  of  the  resultant  Cauchy  and  logarithmic 
singularities  can  be  found  analytically  from  the  basic 
equations  by  assuming  approach  along  the  tangential  path. 
However,  the  strengths  are  in  terms  of  a  distance  along  the 
intersection  from  the  adjacent  vortex  elements  to  the  sheet 
point,  while  the  integrals  in  the  second  terms  of  Eqs.  (2) 
through  (4)  are  over  the  radial  coordinate  along  the  lifting 
line.  This  is  best  handled  by  transforming  the  distance 
along  the  intersection,  as  measured  by  a  curvilinear  co¬ 
ordinate  S  (see  Figure  3),  to  the  radial  coordinate  r0  along 
the  lifting  line.  This  transformation  requires  only  the 
computation  of  the  trajectories  adjacent  to  each  sheet  point. 
The  local  dS/dr0  is  found  by  numerical  differentiation  of 
the  3-point  Lagrangian  interpolation  among  the  3  points. 

The  singularity  strengths  a re  now  completely  determined  in 
terms  of  the  radial  coordinate  r0  along  the  lifting  line. 

The  singulai*  parts  can  now  be  defined  and  are  given  in 
Table  II,  where 

Lu  is  a  function  proportional  to  the  local  curvature 
of  the  projection  of  the  trajectory  in  the  plane 
normal  to  the  x-direction  at  (x,  r,  8), 

Lv  is  the  corresponding  function  for  the  r-direction, 
and 

Lw  is  the  corresponding  function  for  the  0 -direction; 

Q  is  the  component  of  the  input  velocity  normal  to 

the  intersection  of  the  trailing  vortex  sheet  with 
the  plane  through  (x,  r,  8)  normal  to  the  x-axis; 

Cl  is  the  rate  of  change  of  the  input  axial  velocity 
component  u  at  (x,  r,  0)  with  respect  to  the 
parametric  time  t  along  the  trajectory, 

v  is  the  corresponding  derivative  of  v,  and 

w  is  the  corresponding  derivative  of  w;  and 

W  is  the  total  input  velocity  at  the  sheet  point 

(x,  r,  8). 
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SINGULAR  PARTS  OF  INFLUENCE  FUNCTIONS  FOR 
FIELD  POINTS  ALONG  TRAILING  VORTEX  SHEETS 


ruv  +  rvu 


The  singular  parts  of  the  self- induced  velocity  compo¬ 
nents  along  the  trailing  vortex  sheets  can  be  found  by 
analytical  integration  exactly  as  in  INFLOWPT.  Finally, 
the  regular  parts  of  the  velocity  components  are  computed 
exactly  as  in  INFLOWPT.  The  basic  element  spacing  and 
integration  accuracies  were  chosen  the  same  as  with  INFLOWPT 
for  consistency. 


Velocity  Induced  by  Bound  Blade  Vortices 

The  contribution  of  the  bound  blade  vortices  to  the 
induced  velocity  at  field  points  downstream  along  the 
trailing  vortex  sheets  is  in  general  not  zero.  Since  these 
sheet  points  do  not  coincide  with  the  bound  vortex  elements, 
the  appropriate  influence  functions  are  regular  throughout 
the  range  of  interest.  This,  combined  with  the  algebraic 
form  of  these  functions  and  the  fixed  location  of  the 
elements,  makes  the  integration  of  the  first  terms  of  Eqs. 
(2)  through  (4)  relatively  straightforward. 

The  integration  is  carried  out  by  means  of  the  general¬ 
ized  Gaussian  integration  scheme.  The  integrands  of  Eqs. 

(2)  through  (4)  are  computed  at  the  radial  stations  required 
in  the  scheme  by  multiplying  the  calculated  influence  func¬ 
tions  at  these  required  radial  stations  by  the  corresponding 
values  of  r  calculated  from  Eq.  (22).  The  accuracy  used  in 
the  integration  was  -0.00005  of  the  tip  rotational  speed  to 
be  consistent  with  the  velocity  induced  by  the  trailing 
vortex  sheets. 


Checks  of  Computer  Program 

Several  checks  of  SHEETPT  were  made  initially  in  order 
to  gain  confidence  in  the  program  as  well  as  to  begin  inter¬ 
preting  and  using  its  results  to  achieve  a  better  approxi¬ 
mation  to  the  force-free  condition  on  the  trailing  vortex 
sheets.  Unfortunately,  there  were  few  comparable  previous 
results  available  for  checking  purposes.  This  meant  that 
the  checking  and  especially  the  interpretation  of  SHEETPT 
computations  continued  throughout  the  entire  investigation 
and  will  be  treated  subsequently. 

The  only  check  that  could  be  made  was  for  the  special 
case  of  a  propeller  in  cruise  with  the  assumption  of  regular 
helical  trailing  vortex  sheets;  so  we  computed  the  sheet 
velocities  for  the  marine  propeller  that  we  had  considered 
in  the  INFLOWPT  and  MRCHINFO  checks.  Even  here  we  had  no 
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basis  for  comparisons  of  the  total  induced  velocity  on  the 
sheets  although  we  could  check  most  of  the  intermediate 
steps.  For  example,  the  singularity  strengths  computed  from 
the  equations  in  Table  II  reduced  correctly  to  twice  the 
strengths  computed  for  this  propeller  in  INFLOWPT,  according 
to  Eqs.  (24)  through  (26).  The  factor  of  two  arises  because 
the  trajectories  extend  upstream  as  well  as  downstream 
through  the  sheet  points.  More  important,  the  computed 
influence  functions  checked  with  analytical  results  based 
on  the  symmetry  properties  afforded  by  integration  along 
regular  helical  trajectories.  This  check  involved  the  use 
of  MRCHINFO  for  obtaining  integrations  to  finite  distances 
along  the  trajectories.  The  interpretation  and  subsequent 
use  of  the  final  results  for  this  example  will  be  treated  in 
Chapter  4. 

A  number  of  runs  were  then  made  for  a  sample  case  with 
deformation  present.  This  case  was  a  first  approximation 
and  so  was  by  no  means  close  to  a  converged  solution  of  the 
equations.  Nevertheless,  it  was  a  good  example  for  checking 
because  there  was  a  large  amount  of  deformation  present  in 
the  trailing  vortex  sheets.  Most  of  the  procedure  for  cal¬ 
culation  of  the  singular  parts  was  checked  by  hand  since 
these  parts  constitute  the  principal  differences  between 
SHEETPT  and  INFLOWPT.  A  number  of  minor  improvements  were 
made  as  a  result  of  these  detailed  examinations.  We  also 
investigated  the  effects  of  varying  some  of  the  input 
parameters  to  determine  the  best  choices  of  marching  step 
increment  sizes  At  and  locations  of  the  added  trajectories 
that  are  required  for  computing^  and  dS/dr0.  Finally,  in 
order  to  help  explain  certain  features  observed  for  the 
radial  variation  of  the  influence  functions,  a  few  points 
were  recomputed  with  only  the  trailing  vortex  sheet  from  one 
blade  present.  This  gave  a  clear  indication  of  the  strong 
aerodynmaic  interference  among  the  trailing  vortex  sheets, 
especially  downstream,  where  the  deformation  becomes  more 
pronounced. 

From  all  of  these  computations  and  checks,  we  realized 
that  the  complicated  numerical  operations  involved,  coupled 
with  the  large  and  sometimes  locally  severe  deformations  of 
the  sheets,  would  make  careful  investigation  of  all  subse¬ 
quent  results  necessary.  However,  even  if  the  velocity  at 
1  or  2  points  for  each  case  seemed  to  be  unusual  or  difficult 
to  explain,  still,  overall,  the  results  appeared  to  be 
reasonable . 
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CHAPTER  4 


INITIAL  THEORETICAL  RESULTS 


Review  of  Wake  Hypothesis  Development 

The  investigations  described  in  Refs.  10  and  11  were 
begun  with  our  initial  wake  hypothesis  as  an  approximation 
to  the  force -free  condition  of  the  wake.  This  hypothesis 
utilized  the  variations  of  the  three  velocity  components 
with  distance  downstream  along  constant  radii  that  are 
predicted  by  the  theory  of  the  generalized  actuator  disk 
(Ref.  17).  Computations  for  two  different  propellers  indi¬ 
cated  that  the  predicted  performance  was  excessively  opti¬ 
mistic  compared  with  test  data.  Moreover,  preliminary  wake 
velocity  test  data  showed  that  although  the  predicted  con¬ 
traction  pattern  of  the  envelopes  of  the  elements  of  the 
trailing  vortex  sheets  was  realistic,  the  predicted  effec¬ 
tive  "pitch"  variations  of  the  elements  of  the  trailing 
vortex  sheets  were  not.  We  concluded  that  most  of  our 
attention  should  be  directed  toward  the  determination  of  the 
"pitch"  variations. 

The  resulting  refined  wake  hypothesis  was  based  on  the 
initial  hypothesis,  but  changes  were  made  in  an  attempt  to 
overcome  the  above  difficulties.  This  refined  hypothesis, 
or  force-free  approximation,  was  identical  to  the  initial 
hypothesis  inboard  of  the  0.5  radial  station.  Outboard,  the 
axial  velocity  variations  were  applied  along  streamlines 
instead  of  along  constant  radii;  i.e.,  along  the  trajectories 
of  the  elements  of  the  trailing  vortex  sheets.  The  radial 
velocity  variations  were  fixed  by  assuming  that  both  the 
initial  and  refined  hypotheses  should  give  the  same  contrac¬ 
tion  pattern.  After  an  encouraging  preliminary  trial,  the 
refined  hypothesis  was  used  on  a  third  propeller,  and  good 
agreement  with  performance  test  data  was  obtained  after 
three  successive  approximations.  However,  a  subsequent 
fourth  approximation  suggested  that  less  favorable  agreement 
might  occur  for  a  fully  converged  solution.  By  this  time 
we  had  concluded  that  for  whatever  wake  hypothesis  we 
developed,  it  would  be  necessary  to  check  the  assumed 
velocity  components  along  the  trailing  vortex  sheets  by 
actual  computations  and,  where  feasible,  by  corresponding 
test  data. 

This  conclusion  was  reinforced  by  a  set  of  independent, 
unpublished  computations  for  another  propeller  design.  We 
were  unable  to  achieve  a  convergent  inflow  iteration  using 
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the  refined  wake  hypothesis.  Not  only  was  the  axial  inflow 
becoming  unreasonably  high  over  midblade  and  negative  near 
the  tip,  but  also  the  amount  of  contraction  was  becoming 
excessive.  In  turn,  agreement  with  performance  test  data 
grew  poorer  with  each  successive  approximation. 

Once  program  SHEETPT  had  been  written  and  basically 
checked,  we  turned  to  the  four-bladed,  7-foot-diameter  65 
Activity  Factor  (AF)  propeller  at  a  pitch  setting  of  12.6°, 
as  measured  at  the  0.7  radial  station.  The  blade  character¬ 
istics  of  this  propeller  are  given  in  Appendix  II.  The 
airfoil  section  data  appropriate  to  a  tip  speed  of  800  feet 
per  second  and  the  zeroth  approximations  to  the  axial  inflow 
and  blade  circulation  r  were  determined  on  the  basis  of  Refs. 
6  and  12,  the  same  as  for  the  3 (13168A10P3)  propeller  (see 
pp.  37  and  40  of  Ref.  10).  The  radial  inflow  was  guessed, 
and  the  tangential  inflow  was  found  from  the  generalized 
actuator  disk  result,  -Nr/4rrr.  A  first  approximation  to  the 
inflow  and  r  was  calculated  on  the  basis  of  the  refined  wake 
hypothesis  and  displayed  characteristics  of  previous  unsuc¬ 
cessful  iterations.  So,  instead  of  iterating,  the  first  ap¬ 
proximation  to  the  induced  velocity  components  along  the 
trailing  vortex  sheets  was  computed  for  elements  trailing 
from  the  0.4,  0.7,  0.9,  and  0.975  radial  stations. 

The  computed  variations  of  the  axial  induced  velocity 
u  with  the  distance  x  downstream  increased  considerably  more 
rapidly  than  the  assumed  variations  based  on  the  refinod 
wake  hypothesis.  In  fact,  some  appeared  to  reach  asymptotic 
values  greater  than  the  limit  of  2.0  which  is  inherent  in 
regular  helical  wakes.  In  Figure  4  we  have  plotted  the  x 
variation  of  u  computed  along  the  trajectory  from  the  0.9 
radial  station.  This  and  all  subsequent  variations  are 
normalized  by  the  corresponding  inflow  value.  For  comparison 
we  have  also  plotted  the  variations  of  u  assumed  in  the 
refined  wake  hypothesis  and  those  computed  for  the  marine 
propeller  described  in  the  previous  chapter.  It  must  be 
noted  that  the  variations  plotted  are  the  trailing  contribu¬ 
tions  only.  To  these  must  be  added  the  contribution  from 
the  bound  blade  vortices.  This  contribution  is  usually  much 
less  than  101  of  the  trailing  contribution.  The  bound 
contribution  also  fluctuates  in  sign  as  a  function  of 
azimuthal  location  of  the  sheet  point  in  question  and  decays 
rapidly  with  x.  We  have  neglected  the  bound  contributions 
throughout  the  investigation  because  they  add  only  a  higher 
order  refinement  to  the  basic  variations  and  are  significant 
only  when  a  true,  force- free  solution  is  approached  more 
closely  than  we  have  been  able  to  approach  it. 
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Normalized  Axial  Velocity,  u(x/R_,  0 . 9) /u (0 , 0 . 9) 


Nondimens ional  Axial  Coordinate,  r/Rp 

Figure  4.  Comparison  of  Assumed  and  Computed  x 
Variation  of  Axial  Induced  Velocity 
Along  Trajectory  From  r0/Rp  =  0.9, 

65  AF  Propeller,  P0.7R  =  12.6°, 

Refined  Wake  Hypothesis. 
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The  computed  x-variations  of  the  radial  velocity  v 
indicated  initially  that  more  contraction  is  required,  but 
downstream  the  results  were  not  clear  at  all.  The  computed 
x-variations  of  the  tangential  velocity  w  increased  very 
rapidly  to  approach  apparent  asymptotic  values,  which  are 
greater  than  the  mean  value  found  by  conservation  of  the 
circulation.  It  should  be  recalled  that  w  is  small  compared 
to  the  local  propeller  rotational  velocity  nr  and  so  is  of 
secondary  importance. 

On  the  basis  of  the  evidence  we  gathered  through  all 
applications  of  the  refined  wake  hypothesis  as  well  as 
through  the  comparisons  of  the  assumed  and  computed  velocity 
variations,  we  concluded  that  the  refined  wake  hypothesis 
did  not  provide  a  representation  of  the  wake  that  is  general¬ 
ly  adequate.  We  continued  to  feel  strongly  that  the  effec¬ 
tive  "pitch"  of  the  trailing  vortex  system  was  the  most 
important  feature  of  the  wake,  provided  that  a  reasonable 
contraction  pattern  existed.  Moreover,  it  appeared  that  it 
would  be  advisable  to  compute  the  x-variations  of  the  induced 
velocity  components  early  in  the  iteration  to  aid  in  obtain¬ 
ing  an  adequate  "pitch"  representation  and  so  a  force- free 
solution. 


Formulation  of  New  Wake  Representation 

Simultaneously  with  our  conclusions  that  the  refined 
wake  hypothesis  was  inadequate  as  a  general  approximation  to 
the  force-free  condition  of  the  wake,  new  results  were 
obtained  at  Therm  Advanced  Research,  Inc.,  by  M.  D.  Greenberg 
and  A.  L.  Kaskel  (Ref.  16)  for  an  exact,  force-free  actuator 
disk  model  of  a  propeller  or  rotor  in  static  operation,  in 
and  out  of  ground  effect. 

The  model  in  Ref.  16  considers  the  6-average  or,  equiva¬ 
lently,  the  infinite  blade  number  representation  of  a  propel¬ 
ler  or  rotor  with  a  constant  blade  circulation  distribution 
and  zero  hub  radius.  The  trailing  vortex  system  consists 
of  a  line  vortex  along  the  axis  of  rotation  and  a  deformed 
vortex  tube  from  the  tip.  The  exact  kinematic  and  dynamic 
equations  have  been  solved  numerically  for  the  entire  flow 
field.  This  formulation  has  the  property  that,  if  the 
tangential  velocity  w  is  neglected,  the  entire  streamline 
contraction  pattern  is  analytically  independent  of  the  total 
loading.  If  w  is  included,  the  pattern  is  numerically 
virtually  independent  of  the  total  loading.  This  contraction 
pattern  out  of  the  ground  effect,  illustrated  in  Figure  5, 
agrees  reasonably  well  with  available  flow  visualization 
results . 
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Figure  5.  Contraction  Pattern  of  Envelopes  of  Trailing  Vortex  Elements 
New  Wake  Representation  Based  on  Ref.  16. 


The  first  important  feature  of  our  new  wake  representa¬ 
tion  is  the  specification  of  the  contraction  pattern  by  these 
actuator  disk  results;  i.e.,  the  ratio  of  the  radial- to-axial 
velocity  components  v/u  is  fixed  along  each  element  of  the 
trailing  vortex  sheets  by  Ref.  16.  This  is  incorrect  far 
inboard  because  the  radial  velocity  is  usually  positive  there 
due  to  the  effect  of  the  distributed  trailing  vorticity,  but 
this  region  has  very  little  effect  outboard.  Also,  the  x- 
variations  of  the  tangential  velocity  w  along  all  the  trajec¬ 
tories  are  fixed  according  to  the  values  computed  for  the 
65  AF  propeller  at  the  12.6°  pitch  setting.  This  leaves  the 
x-variations  of  the  axial  velocity  u,  which  determine  the 
effective  "pitch"  along  the  trajectories  of  the  elements  of 
the  trailing  vortex  sheets,  as  the  principal  variables  at 
our  disposal  for  approximating  a  force-free  wake. 

This  representation  of  the  wake  insured  in  advance  a 
reasonable  contraction  pattern  and  allowed  us  to  concentrate 
on  working  toward  the  force-free  effective  "pitch."  However, 
v  is  computed  along  with  the  other  components,  and  this 
permitted  us  to  draw  conclusions  afterwards  about  the  valid¬ 
ity  of  the  assumed  fixed  contraction  pattern. 


65  AF  Propeller,  12.6°  Pitch  Setting  -  Initial  Computations 

The  first  case  treated  by  means  of  the  new  wake  repre¬ 
sentation  was  again  the  12.6°  pitch  setting  of  the  65  AF 
propeller.  We  started  from  the  same  zeroth  approximations 
and  based  the  x-variations  of  the  axial  velocity  u  on  the 
computations  for  this  case  with  the  refined  wake  hypothesis. 

Three  successive  approximations  to  the  inflow  and  blade 
circulation  were  computed,  and  the  axial  inflow,  although 
initially  reasonable,  was  beginning  to  diverge  as  we  had 
observed  for  some  of  the  refined  wake  hypothesis  computations. 
Therefore,  we  computed  the  third  approximation  to  the  x- 
variations  of  the  velocity  components  along  the  trailing 
vortex  sheets.  The  x-variations  of  u  increased  less  rapidly 
with  x  than  assumed  inboard  of  the  0.9  radial  station,  but 
at  the  0.975  station  the  x-variations  increased  much  more 
rapidly  and  reached  a  higher  asymptotic  value. 

The  outboard  variations  were  modified  accordingly  to 
agree  with  the  computed  values,  and  the  inflow  was  recomputed 
with  the  second  approximation  as  input.  The  resulting  axial 
inflow  seemed  to  be  too  low  over  midblade  and  too  high  at 
the  tip.  We  decided  that  when  modifying  the  variations,  we 
should  not  make  changes  to  the  full  extent  indicated  by  the 
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computed  variations.  Therefore,  we  changed  the  variations 
to  lie  roughly  halfway  between  the  previous  two  sets  and 
again  recomputed  the  inflow  with  the  second  approximation 
as  input.  The  axial  inflow  seemed  to  be  reasonable,  but 
the  recomputed  x-variations  of  u  were  close  to  those  assumed 
initially.  We  repeated  this  cycle  of  changing  the  x-varia¬ 
tions  of  u  and  recomputing  the  inflow  and  x-variations  once 
more.  A  convergent  trend  did  seem  to  be  emerging  in  the  x- 
variations  of  u,  but  the  inflow  and  blade  circulation  were 
still  not  close  to  convergence.  Furthermore,  the  entire 
iterative  process  was  very  inefficient  and  unwieldy,  so  we 
suspended  calculations  for  this  case.  The  calculations  were 
resumed  later  and  are  described  in  Chapter  6. 


65  AF  Propeller,  8.2°  Pitch  Setting 

We  turned  our  attention  next  to  the  8.2°  pitch  setting 
of  the  same  65  AF  propeller  operating  at  a  tip  speed  of  800 
feet  per  second.  This  case  is  significant  because  the  in¬ 
stantaneous  velocity  components  in  the  propeller  wake  had 
been  measured  by  means  of  hot-wire  anemometry  in  the  experi¬ 
mental  portion  of  this  investigation.  These  data  were 
obtained  to  guide  and  check  the  theoretical  computations. 

The  azimuthal  variations  of  the  axial  and  radial  velocity 
components  at  several  points  in  the  wake  are  presented  in 
Appendix  III,  while  the  full  details  of  the  acquisition  and 
reduction  of  the  data  are  contained  in  the  test  report 
prepared  by  D.  C.  Gilmore  and  I.  S.  Gartshore  of  McGill 
University,  Ref.  14. 

The  theoretical  calculations  began  by  obtaining  zeroth 
approximations  to  the  inflow  and  blade  circulation  in  the 
same  way  as  for  the  12.6°  pitch  setting  of  this  propeller. 

We  selected  the  same  x-variations  of  u  along  the  trailing 
vortex  sheets  that  we  had  used  originally  for  the  12.6° 
setting  and  computed  a  first  approximation  to  the  inflow. 

The  axial  inflow  was  nearly  50%  higher  over  midblade 
but  was  strongly  negative  outboard  of  the  0.95  radial  station. 
This  held  no  promise  for  a  successful  iteration,  so  we 
modified  the  x-variations  of  u  outboard  of  the  0.95  station 
to  increase  more  rapidly  with  x  and  to  reach  a  higher  as¬ 
ymptotic  value.  The  inflow  was  then  recomputed  but  still 
held  no  promise  for  iteration  because  the  axial  inflow  was 
essentially  unchanged  outboard. 

We  examined  the  wake  velocity  data  next  to  assess  their 
use  as  a  guide  for  achieving  convergent  theoretical  compu¬ 
tations.  A  principal  feature  of  the  measurements  was  that 
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a  hot-wire  signal  which  repeated  with  each  blade  passage 
(see  Figure  6,  for  example}  could  be  obtained  only  in  a 
limited  region  behind  the  middle  of  the  blade.  On  the  other 
hand,  inboard  near  the  hub  and  outboard  near  the  tip,  the 
signal  was  apparently  completely  random  (see  Figure  7)  and 
did  not  repeat  for  each  blade  passage.  The  extent  and  growth 
downstream  of  the  nonrepeatable  region  are  shown  in  Figure  8 
for  this  case,  as  well  as  for  the  1 20  AF  propeller  at  the 
10.0°  setting;  see  Chapter  6.  Beyond  the  outer  limit  of  this 
region  and  upstream  of  the  propeller,  a  periodic  signal  was 
again  obtained.  We  have  been  unable  to  explain  the  cause  of 
the  random  signals. 

The  randomness  precluded  a  meaningful  resolution  of  the 
basic  data  into  the  instantaneous  velocity  components,  and 
so,  practically,  we  do  not  have  useful  data  in  the  tip  region. 
In  our  theoretical  computations,  it  is  the  "pitch”  of  the 
trailing  vortex  sheets  in  this  tip  region  that  is  most  impor¬ 
tant,  because  the  tip  elements  principally  control  the  level 
of  the  inflow  over  the  rest  of  the  blade.  Thus,  the  absence 
of  data  in  the  tip  region  limits  the  use  of  the  wake  data 
to  checking  the  computed  values  over  midblade;  that  is,  the 
data  can  be  used  as  a  check  but  not  as  a  guide  in  defining 
the  "pitch"  where  the  "pitch"  is  most  important.  Therefore, 
in  this  case  where  we  were  unable  to  iterate  successfully, 
the  data  were  not  of  immediate  use. 

The  principal  difficulty  in  iterating  appears  to  be  an 
underloading  of  the  propeller  near  the  tip.  This  was  ob¬ 
served  in  the  zeroth  approximation  where  the  blade  circu¬ 
lation  r  was  very  low.  The  indicated  negative  axial  inflow 
near  the  tip  in  the  first  approximation  would  lead  to  a 
larger  r  there,  but  we  were  unable  to  treat  negative  axial 
velocities  with  any  confidence.  When  the  axial  velocities 
are  very  small,  as  in  the  zeroth  approximation,  or  negative, 
as  in  the  first,  the  exact  location  of  the  trailing  vortex 
elements  is  extremely  important  because  small  changes  in 
their  position  can  cause  large  interference  effects,  as  we 
found  in  earlier  studies  (see  pp.  31  and  33  of  Ref.  10). 

In  view  of  the  difficulties  encountered  in  our  theo¬ 
retical  computations,  we  changed  our  approach  significantly 
at  this  point  in  order  to  overcome  these  difficulties.  The 
computations  toward  this  end  are  described  in  the  next 
chapter . 
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CHAPTER  5 


THEORETICAL  RESULTS  FOR  REPRESENTATIVE  BLADE  LOADINGS 


General 


An  approach  that  an  experienced  propeller  designer 
might  find  useful  as  an  aid  in  th  development  of  improved 
blade  designs  can  be  outlined  as  follows.  First,  a  desired 
blade  circulation  distribution  r  or,  equivalently,  a  blade 
loading  distribution  would  be  assumed.  The  solution  for 
the  corresponding  inflow  components  would  then  be  found  and 
used  in  conjunction  with  r  to  determine  an  actual  blade 
design.  The  utility  of  this  approach  rests  ultimately,  of 
course,  on  a  proven,  suitable  method  for  calculating  the 
inflow;  i.e.,  a  method  that  is  based  on  a  force-free  wake 
and  that  has  been  verified  by  comparisons  with  test  data  for 
widely  varying  existing  propeller  designs. 

In  the  stage  of  our  theoretical  development  when 
computations  for  representative  blade  loadings  were  under 
taken, such  computations  could  only  serve  to  demonstrate  the 
feasibility  of  such  an  approach.  On  the  other  hand,  the  same 
approach  offered  an  important  simplification  computationally, 
since  r  remains  fixed  throughout  the  iteration  and  Eq.  (13) 
is  not  needed  in  conjunction  with  Eqs.  (2)  through  (4)  to 
find  the  inflow.  We  pursued  this  approach  at  this  point  of 
the  investigation  because  of  the  promise  it  offered  in  pro¬ 
viding  important  experience  and  insight  into  the  attainment 
of  force-free  solutions  to  Eqs.  (2)  through  (4). 

We  selected  both  a  four-bladed  r,  called  Representative 
Blade  Loading  One  (RBL  1),  and  a  three-bladed  r,  called  RBL 
2.  Each  of  these  distributions  was  found  by  extrapolating 
T.  Theodorsen's  moderate  loading  extension  of  the  Goldstein 
optimum  analysis  (see  Eq.  (1)  on  p.  23  of  Ref.  29)  to  zero 
advance.  Contraction  was  ignored  in  this  extrapolation,  and 
each  r  was  truncated  to  a  hub  radius  at  the  O.J  radial 
station.  Both  distributions  were  selected  to  have  comparable 
loadings  and  are  based  upon  a  "displacement  velocity  at 
infinity  downstream"  (see  Ref.  29)  of  1/6  of  the  propeller 
tip  speed.  The  corresponding  induced  thrust  coefficients 
(found  by  eliminating  the  lift  coefficient  between  Eqs.  (13) 
and  (14)  while  simultaneously  neglecting  profile  drag  losses; 
i.e.,  y  =  0)  are  0.161  for  RBL  1  and  0.154  for  RBL  2.  The 
circulation  distributions  for  RBL  1  and  RBL  2  are  plotted 
later  along  with  the  computed  inflow  in  Figure  9  and  11, 
respectively. 
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The  zeroth  approximations  to  the  axial  and  tangential 
components  of  the  inflow  were  found  from  the  equations 
implicit  in  the  analysis  of  M.  Knight  and  R.  Hefner  (Ref. 
19).  In  particular,  the  axial  inflow  is  the  square  root  of 
Nnr/4ir,  and  the  tangential  inflow  is  the  familiar  -Nr/4rr. 
However,  we  faired  both  components  to  finite  values  at  the 
tip  to  be  more  realistic. 


Representative  Blade  Loading  One 

The  four-bladed  case,  RBL  1,  was  considered  first.  With 
the  blade  circulation  r  fixed  and  the  zeroth  approximations 
to  the  inflow  determined  as  above,  it  remained  only  to  speci¬ 
fy  the  x-variations  of  the  axial  velocity  u  along  the  trajec¬ 
tories.  We  chose  variations  which  increased  with  x  initial¬ 
ly  off  the  blade  at  rates  consistent  with  our  preceding  com¬ 
putations  but  which  reached  asymptotic  values  of  exactly 
2.0.  These  variations,  then,  are  not  greatly  different  from 
those  assumed  in  the  refined  wake  hypothesis  and  so  have  some 
basis  in  our  previous  results. 

Four  successive  approximations  to  the  i* f low  were 
computed.  The  axial  iteration  factors  used  with  Eq.  (19) 
differed  from  one  approximation  to  the  next  as  we  proceeded 
cautiously  toward  a  solution.  We  followed  this  procedure 
throughout  the  computations  for  RBL  1  and  RBL  2.  Tangential 
iteration  factors  of  1.0  were  generally  taken  throughout. 

The  fourth  approximation  to  the  axial  inflow,  with  an  itera¬ 
tion  factor  of  1.0,  converged  to  within  ±4.0$  of  the  third 
approximation  everywhere  along  the  blade  and  to  within  ±1.1$ 
outboard  of  the  0.9  radial  station.  We  judged  this  to  be 
sufficient  convergence  for  checking  the  x-variations  of  u 
downstream.  The  fourth  approximations  to  the  axial  and 
tangential  inflow  are  plotted  in  Figure  9  along  with  r  and 
the  zeroth  approximations. 

In  Figure  10  are  plotted  the  assumed  radial  inflow  used 
as  an  input  to  the  calculation  of  the  fourth  approximation; 
i.e.,  as  fixed  by  the  assumed  contraction  pattern  and  the 
third  approximation  to  the  axial  inflow,  and  the  computed 
fourth  approximation  with  a  radial  iteration  factor  of  1.0. 

The  contraction  pattern  is  clearly  inadequate  inboard,  as 
we  knew  in  advance.  However,  the  agreement  outboard  is 
fairly  good  except  that  near  the  tip  the  indications  are  that 
the  fixed  pattern  is  providing  too  much  contraction.  Finally, 
the  computed  induced  power  coefficient  (found  by  eliminating 
the  lift  coefficient  between  Eqs.  (13)  and  (15)  with  y  =  0) 
is  1.0476  and  leads  to  an  induced  figure  of  merit  of  1081. 
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Figure  9.  Blade  Circulation  and  Axial  and  Tangential 
Inflow,  Representative  Blade  Loading  One, 
First,  Second,  and  Third  Assumed  x- Variations 
of  Axial  Velocity. 
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The  next  step  was  to  compute  the  fourth  approximation 
to  the  x-variations  of  the  induced  velocity  components  along 
the  trailing  vortex  sheets.  At  the  0.4  and  0.7  radial  sta¬ 
tions  the  computed  x-variations  of  u  increased  less  rapidly 
with  x  than  assumed,  while  at  0.9  and  0.975  they  increased 
more  rapidly.  In  light  of  our  experience  with  the  65  AF 
propeller,  we  changed  the  variations  by  only  about  one-third 
of  the  differences  between  the  assumed  and  computed  varia¬ 
tions.  This  was  performed,  moreover,  in  two  steps.  First, 
we  changed  the  variations  outboard  of  the  0.8  station  only 
and  computed  the  inflow.  Then  we  changed  the  variations  in¬ 
board  of  0.8  as  well  and  recomputed  the  inflow.  The  axial 
inflow  between  the  two  differed  by  less  than  1.0%  over  the 
entire  blade.  This  gave  valuable  confirmation  of  our  earlier 
conclusions  that  it  is  the  outboard  variations  that  need  to 
be  represented  more  accurately. 

With  these  x-variations  as  finally  modified,  and  begin¬ 
ning  with  the  previous  third  approximation  as  the  zeroth, 
two  successive  approximations  to  the  inflow  were  computed. 

The  axial  inflow,  which  converged  to  within  -1.7%  everywhere 
and  within  -0.2%  outboard  of  the  0.9  station,  is  also  plotted 
in  Figure  9.  The  assumed  and  computed  radial  inflow  compari¬ 
son  was  somewhat  poorer  than  for  the  first  set  of  variations. 
The  induced  performance  became  even  less  reasonable,  too, 
with  a  power  coefficient  of  0.0434  and  a  figure  of  merit  of 
119%. 

The  corresponding  second  approximation  to  the  x-varia¬ 
tions  of  the  induced  velocity  components  was  computed.  Sur¬ 
prisingly,  the  computed  x-variations  of  u  were  virtually 
identical  to  those  computed  with  the  first  set  of  assumed 
variations,  so  we  modified  them  to  the  full  extent  indicated. 

Assuming  this  third  set  of  x-variations  of  u  and  taking 
the  immediately  previous  first  approximation  as  the  zeroth, 
three  successive  approximations  to  the  inflow  were  computed. 
The  axial  inflow,  again  shown  in  Figure  9,  converged  to 
within  -1.5%  everywhere  along  the  blade.  The  comparison  of 
the  assumed  and  computed  radial  inflow  is  presented  in 
Figure  10,  and  agreement  is  even  poorer  than  for  the  first 
two  sets  of  x-variations  of  u.  Furthermore,  the  indications 
are  clearer  that  the  fixed  pattern  provides  too  much  contrac¬ 
tion.  The  trend  to  less  reasonable  induced  performance 
continued  with  a  power  coefficient  of  0.0355  and  a  figure 
of  merit  of  1451. 

Both  the  second  and  third  approximations  to  the  x- 
variations  of  the  induced  velocity  components  were  computed 
with  this  assumed  set  of  variations.  The  axial  inflow 


differed  between  these  two  approximations  by  less  than  ^5.0% 
and  then  only  in  the  sense  that  the  second  did  not  have  the 
minimum  at  the  0.9  radial  station  (See  Figure  9)  but  Was 
instead  nearly  constant  outboard  of  0.8.  At  the  0.4  and  0.7 
radial  stations,  the  computed  x- variations  of  u  were  nearly 
identical  in  the  two  cases  and  were  very  close  to  the  as¬ 
sumed  variations.  However,  at  the  0.9  station  especially, 
and  to  a  lesser  degree  at  the  0.975  station,  the  computed 
x-variations  of  u  were  markedly  different  in  the  two  approxi¬ 
mations.  Both  approximations  increased  with  x  at  close  to' the 
assumed  rates  initially .  However ,  the  second  approximations  con¬ 
tinued  to  increase  at  close  to  the  assumed  rates, while  the  third 
approximations  reached  maxima  at  about  0.2  blade  radius  down¬ 
stream  and  then  decreased.  The  apparent  cause  of  these 
large  differences  is  the  pronounced  change  in  the  degree  of 
distortion  of  the  trailing  vortex  sheets,  resulting  from 
differences  in  axial  inflow  which  are  relatively  small  in 
magnitude  but  which  vary  rapidly  with  radius.  In  particular, 
these  gradients  in  axial  velocity  outboard  led  to  a  local 
stretching  out  of  the  trailing  vortex  sheets  in  the  third 
approximation  that  was  not  present  in  the  second.  Therefore, 
for  conservation  of  circulation,  the  trailing  vorticity  out¬ 
board  is  more  thinly  distributed  over  the  sheet  in  the  third 
approximation.  Consequently,  the  induced  velocity  is  con¬ 
siderably  smaller. 

These  computations  completed  consideration  of  RBL  1. 

We  will  discuss  the  results  more  fully  after  a  description 
of  the  RBL  2  computations  which  were  carried  out  next. 


Representative  Blade  Loading  Two 


Similar  computations  were  made  for  RBL  2  with  the  blade 
circulation  and  the  zeroth  approximations  to  the  inflow  based 
on  Ref.  19,  as  described  earlier.  To  provide  comparisons 
with  RBL  1,  we  selected  the  same  set  of  x-variations  of  the 
axial  velocity  u  that  we  had  used  initially  for  RBL  1. 


Three  successive  approximations  to  the  inflow  were 
computed.  The  convergent  trend  was  rapid,  and  the  third 
approximation  to  the  inflow  was  within  -1.01  of  the  second. 
The  axial  and  tangential  inflows  are  plotted  in  Figure  11 
along  with  r  and  the  zeroth  approximations.  Qualitatively, 
the  assumed  and  computed  radial  inflows  differed  like  they 
did  in  the  third  set  of  x-variations  of  u  for  RBL  1.  The 
induced  performance  is  even  less  realistic  than  the  compa¬ 
rable  RBL  1  case,  with  a  power  coefficient  of  0.0406  and  a 
figure  of  merit  of  1191. 
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Figure  11.  Blade  Circulation  and  Axial  and  Tangential 
Inflow,  Representative  Blade  Loading  Two, 
First  and  Second  Assumed  x-Variations  of 
Axial  Velocity. 


Next,  the  third  approximation  to  the  x-variations  of 
the  induced  velocity  components  along  the  trailing  vortex 
sheets  was  computed.  At  the  0.4  and  0.7  radial  stations, 
the  computed  x-variations  of  7  increased  Jess  rapidly  with 
x  than  assumed,  just  as  for  RBL  1.  At  0.9,  the  computed 
variation  increased  more  rapidly  than  assumed  initially  and 
then  leveled  off  at  a  higher  asymptotic  value.  At  0.975, 
the  computed  variation  also  increased  more  rapidly  initially 
but  then  reached  a  maximum  and  subsequently  decreased;  see 
Figure  12.  This  is  v.he  same  behavior  that  we  observed  for 
RBL  1.  We  then  modified  the  x-variations  of  u  to  be  identi¬ 
cal  with  the  tlird  set  of  variations  for  RBL  1  inboard  of 
the  0.85  station  and  to  follow  the  newly  computed  RBL  2 
variations  outboard  of  0.85,  simply  ignoring  the  decrease 
after  the  maximum;  see  Figure  12. 

With  the  second  set  of  x-variations  of  u,  and  with  the 
previous  second  approximation  as  the  zeroth,  two  successive 
approximations  to  the  inflow  were  computed.  The  second  ap¬ 
proximation  to  the  axial  inflow,  which  converged  to  within 
+2.6$  of  the  first  approximation  inboard  of  the  0.7  radial 
station  and  within  +0.1$  outboard,  is  also  plotted  in  Figure 
11.  The  assumed  and  computed  radial  inflow  comparison  was 
poorer  yet,  continuing  the  deteriorating  trend  that  was  also 
shown  by  the  RBL  1  cases  in  Figure  10.  The  induced  perfor¬ 
mance  became  less  reasonable,  too,  with  a  power  coefficient 
of  0.0357  and  a  figure  of  merit  of  135$. 

The  corresponding  second  approximation  to  the  x-varia¬ 
tions  of  the  induced  velocity  components  was  computed.  At 
the  0.4  and  0.7  stations,  the  computed  x-variations  of  u 
were  in  good  agreement  with  the  assumed  values.  At  0.9,  the 
computed  variation  agreed  well  initially  with  that  assumed 
but  reached  a  higher  asymptotic  value.  At  0.975,  the  x- 
variation  of  u  also  agreed  well  initially  but  then  increased 
much  more  rapidly  with  x;  see  Figure  12.  It  is  highly 
encouraging  that  a  maximum  did  not  appear,  at  least  as  far 
downstream  as  computations  were  made.  This  appears  to  indi¬ 
cate  that  the  maxima  and  subsequent  decreases  are  only 
temporary  phenomena  which  may  occur  at  particular  steps 
toward  a  force-free  solution  due  to  local  peculiarities  in 
the  trailing  vortex  sheet  distortion.  The  peculiarity  here 
again  seems  to  be  a  result  of  the  relative  stretching  of  the 
sheets,  which  is  considerably  less  than  it  was  for  the  ini¬ 
tial  set  of  variations. 

It  is  clear  from  the  results  presented  for  RBL  1  and 
RBL  2  that  we  have  not  yet  obtained  force-free  solutions  for 
these  blade  loadings.  Nevertheless,  the  successful  inflow 
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iterations  that  we  have  achieved  led  us  to  return  to  the 
specific  propeller  designs  so  that  comparisons  could  be 
made  with  test  data  wherever  possible.  Before  these  are 
reported,  though,  a  number  of  useful  observations  can  be 
made  from  the  RBL  results. 


Observations  for  Representative  Blade  Loadings 

We  are  presenting  in  this  section  certain  observations 
based  upon  our  computations  for  RBL  1  and  RBL  2.  However, 
it  must  be  recognized  that  these  loadings  are  indicative  of 
thrust  coefficients  corresponding  to  the  higher  blade  pitch 
settings  of  typical  propeller  designs.  Thus,  the  applica¬ 
bility  of  these  observations  and  the  subsequent  generaliza¬ 
tions  that  we  draw  are  unexplored  for  loadings  at  thrust 
coefficients  more  typical  of  lower  pitch  settings.  Regard¬ 
less,  we  o^fer  them  in  the  hope  that  they  might  provide  aid 
in  initiating  computations  for  such  cases. 

The  first  area  of  observations  concerns  an  approximate 
correlation  that  we  have  found  for  the  axial  inflow.  With 
the  blade  circulation  distribution  r  fixed,  the  singular 
part  of  the  axial  inflow  is  essentially  fixed,  too,  because 
the  coefficients  of  the  singular  parts  in  Eq.  (24)  are 
virtually  independent  of  the  inflow,  as  we  have  observed 
before  (see  p.  31  of  Ref.  10).  As  an  iteration  proceeds, 
then,  for  a  given  set  of  x-variations  of  the  axial  velocity 
u,  it  is  only  the  regular  part  of  the  axial  inflow  that  dif¬ 
fers  from  one  approximation  to  the  next.  By  consideration 
of  each  application  of  the  Biot-Savart  operator  in  Eq.  (2), 
we  found  that  the  product  of  the  input  axial  inflow  at  a 
representative  radial  station  near  the  tip,  0.975  in  particu¬ 
lar,  times  the  regular  part  of  the  output  axial  inflow  at 
any  radial  station  was  relatively  invariant  throughout  the 
iteration.  In  an  approximate  sense,  this  correlation  can  be 
explained  by  consideration  of  the  Bio-Savart  operator  in  Eq. 
(2)  as  broken  into  singular  and  regular  parts.  The  corre¬ 
lation  seems  to  imply  that  the  regular  part  of  the  trailing 
axial  influence  function  at  any  radial  station  is  approxi¬ 
mately  a  function  of  radius  only  divided  by  the  input  axial 
inflow  at  the  0.975  station,  since  multiplication  of  both 
sides  of  the  equation  by  this  quantity  gives  an  integral 
which  is  essentially  invariant  from  one  approximation  to  the 
next.  The  validity  of  the  inverse  proportionality  is  sug¬ 
gested  by  the  generalized  actuator  disk  result  (see  Eq.  22 
of  Ref.  17),  that  the  mean  axial  inflow  at  any  radius  is 
inversely  proportional  to  the  free-stream  velocity,  which 
determines  the  pitch  of  the  trailing  vortex  sheets  in  that 
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case.  In  the  present  static  case,  the  input  axial  inflow 
at  0.975  is  a  measure  of  the  effective  "pitch”  of  the  trail¬ 
ing  vortex  sheets.  Therefore,  the  correlation  is  a  quanti¬ 
tative  verification  of  our  persistent  observation  that  the 
effective  "pitch"  of  the  trailing  vortex  sheets  near  the  tip 
controls  the  magnitude  of  the  axial  inflow  at  all  radii. 

The  second  area  of  observations  concerns  the  inflow 
iterations  for  a  specified  r  and  set  of  x-variations  of  u. 
These  iterative  properties  have  emerged  clearly  from  all  the 
RBL  1  and  RBL  2  computations.  For  both  loadings,  the  zeroth 
approximations  to  the  axial  and  tangential  inflow  based  on 
Ref.  19  were  adequate.  Axial  inflow  convergence  always 
began  outboard  and  then  continued  to  improve  farther  in¬ 
board  with  each  successive  approximation.  This  results  from 
the  general  property  that  the  velocity  induced  by  a  distorted 
helical  vortex  is  much  larger  at  radii  smaller  than  itself 
than  at  radii  larger  than  itself.  The  tangential  inflow 
converged  rapidly  but  is  different  outboard  from  -Nr/4Trr, 
as  can  be  seen  from  Figures  9  and  11. 

Axial  inflow  convergence  proceeds  very  rapidly  when  the 
x-variations  of  u  are  changed.  In  particular,  the  axial 
inflow  outboard  of  the  0.9  radial  station  changes  only 
slightly  (see  Figures  9  and  11),  although  the  inboard  sta¬ 
tions  are  affected  greatly.  This  demonstrates  again  the 
insensitivity  of  the  axial  inflow  to  local  changes  in  the 
x-variation  but  the  profound  sensitivity  inboard  to  x- 
variation  changes  outboard. 

On  the  basis  of  these  observations,  we  can  recommend  a 
quick  way  to  achieve  axial  inflow  convergence  for  a  fixed  r. 
First,  compute  the  Biot-Savart  operators  with  the  zeroth 
approximation  as  input.  Use  an  axial  iteration  factor  of 
1.0  to  determine  the  first  approximation  outboard  of  the  0.8 
radial  station.  With  the  resulting  axial  inflow  at  0.975, 
use  the  correlation  outlined  above  to  determine  the  axial 
inflow  at  the  remaining  radii.  For  all  subsequent  approxi¬ 
mations,  use  axial  iteration  factors  of  1.0  over  the  entire 
blade.  Although  the  above  procedure  might  minimize  computer 
time,  overall  it  might  be  just  as  efficient  to  use  axial 
iteration  factors  of  1.0  over  the  entire  blade  for  the  first 
approximation  as  well  as  the  subsequent  ones.  Moreover,  when 
the  x-variations  of  u  are  changed,  axial  iteration  factors 
of  1.0  should  be  used  because  the  axial  inflow  changes  only 
slightly  near  the  tip  and  so  converges  very  quickly  elsewhere. 

The  third  area  of  observations  concerns  the  iteration 
toward  a  satisfactory  force-free  approximation  by  means  of 
changes  in  the  x-variations  of  u.  Both  RBL  1  and  RBL  2  have 
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similar  properties  in  this  regard.  There  seems  to  be  a 
definite  trend  toward  agreement  between  the  assume!  and 
computed  variations  immediately  downstream  of  the  lifting 
line.  The  results  in  Figure  12  are  typical  of  this  trend. 
Farther  downstream,  however,  the  computed  x-variations  of  u 
may  have  features  like  the  observed  maxima,  which  can  result 
even  from  small  changes  in  the  axial  inflow.  There  is  evi¬ 
dence,  as  in  Figure  12,  that  these  are  only  transient 
phenomena  which  disappear  as  we  approach  more  closely  a 
force- free  wake. 

Although  the  x-variations  of  u  appeared  to  be  slowly 
tending  toward  the  force- free  state,  the  computed  radial  in¬ 
flow  was  simultaneously  indicating  poorer  agreement  with  the 
assumed  distribution  given  by  the  specified  contraction  pat¬ 
tern,  as  in  Figure  10.  It  appeared  that  less  contraction 
was  necessary  than  that  assumed.  The  performance  was  also 
becoming  less  realistic  with  improvement  in  the  x-variations 
of  u.  At  the  end  of  our  computations  for  RBL  1  and  RBL  2, 
then,  the  problems  involved  in  iterating  to  a  truly  force- 
free  solution  still  had  to  be  resolved. 


CHAPTER  6 


THEORETICAL  RESULTS  FOR  SPECIFIC  PROPELLERS 


General 


We  returned  to  computations  for  specific  propeller 
designs  at  this  point  after  having  achieved  a  good  deal  of 
experience  through  the  treatment  of  the  two  representative 
blade  loadings.  In  the  present  chapter,  we  shall  discuss 
these  theoretical  results  and  shall  make  comparisons  with 
test  data  whenever  possible  to  aid  in  assessing  the  theo¬ 
retical  development. 

The  first  two  cases  discussed  in  this  chapter  are 
examples  of  high  pitch  settings  and  high  thrust  coefficients. 
Consequently,  they  have  blade  circulation  and  inflow  distri¬ 
butions  which  resemble  the  RBL  most  closely.  These  cases 
could  be  iterated  successfully.  After  they  had  been  brought 
to  a  reasonable  degree  of  completion,  including  a  change  in 
the  x-variations  of  u  for  one  of  them,  we  pursued  calcula¬ 
tions  at  lower  pitch  settings.  Meaningful  results  were 
obtained  for  two  of  these,  which  are  indicative  of  an  impor¬ 
tant  regime  of  static  propeller  operation. 

At  all  times  during  these  computations,  we  fed  in  our 
best  previous  experience,  even  if  we  had  no  assurance  that 
a  particular  choice  (say,  of  iteration  factors,  apparent  cor¬ 
relations,  etc)  was  a  valid  generalization  of  sn  earlier 
result.  Also,  although  the  cases  are  described  basically  in 
chronological  order,  some  of  them  were  carried  out,  at  least 
in  part,  simultaneously. 


120  AF  Propeller,  17.0°  Pitch  Setting 

The  first  case  considered  after  the  RBL  was  the  17.0° 
pitch  setting,  as  measured  at  the  0.7  radial  station,  of  the 
four-bladed,  7-foot-diameter,  120  Activity  Factor  (AF)  pro¬ 
peller.  The  characteristics  of  this  blade  design  are  given 
in  Appendix  II.  This  design  is  simply  a  scaled-up  version 
of  the  basic  65  AF  blade,  although  it  has  a  different  blade 
twist  distribution  8.  The  airfoil  section  data  appropriate 
to  a  tip  speed  of  700  feet  per  second  and  the  zeroth  approxi¬ 
mation  to  the  axial  inflow  and  blade  circulation  distribution 
r  were  determined  on  the  basis  of  Refs.  6  and  12  exactly  as 
for  the  65  AF  propeller.  The  zeroth  approximation  to  the 
tangential  inflow  was  found  as  usual  from  r  and  the  general- 
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ized  actuator  disk  result  of  -Nr/4irr.  The  x-variations 
downstream  of  the  axial  velocity  u  were  chosen  to  be  the 
same  set  that  was  used  initially  for  RBL  1  and  RBL  2. 

Seven  successive  approximations  to  r  and  the  inflow 
were  computed  to  achieve  convergence  of  the  seventh  approxi¬ 
mation  to  the  axial  inflow  to  within  ±1.6%  of  the  sixth 
approximation.  The  axial  iteration  factors  used  differed 
from  one  approximation  to  the  next  as  we  tried  to  proceed 
quickly  to  convergence.  The  seventh  approximations  to  r 
and  the  axial  and  tangential  inflow  are  plotted  in  Figure  13 
along  with  the  zeroth  approximations.  The  assumed  and 
computed  radial  inflow  compared  favorably  here,  better  than 
in  any  of  the  RBL  cases.  However,  the  predicted  performance 
is  too  optimistic  compared  with  the  data  acquired  at  this 
pitch  setting  in  our  test  program;  see  Table  III. 

Next,  the  seventh  approximation  to  the  x-variations  of 
the  induced  velocity  components  along  the  trailing  vortex 
sheets  was  computed.  At  the  0.4  and  0.7  radial  stations, 
the  computed  x-variations  of  u  increase  less  rapidly  than 
assumed,  comparable  to  RPL  1  and  RBL  2.  At  0.9,  the  computed 
variation  agrees  closel>  with  the  assumed  variation.  At 
0.975,  the  computed  variation  increases  much  more  rapidly 
than  assumed  and  has  a  maximum  and  subsequent  decrease. 

The  maximum  is  higher  than  for  either  RBL  1  or  RBL  2.  We 
modified  the  x-variations  of  u  accordingly,  but  we  ignored 
the  decrease  beyond  the  maximum  at  the  0.975  radial  station. 

Four  successive  approximations  to  r  and  the  inflow  were 
computed  with  these  modified  x-variations  of  u  and  with  the 
previous  sixth  approximation  as  the  zeroth.  The  axial  inflow 
converged  to  within  ±1.8'<  except  at  the  0.2  radial  station, 
where  blade  stall  apparently  prevented  convergence.  Fewer 
approximations  were  necessary  here  because  r  and  the  axial 
inflow  outboard  did  not  change  greatly;  see  Figure  13.  The 
assumed  and  computed  radial  inflow  comparison  was  appreciably 
less  favorable  than  it  had  been  initially,  with  indications 
that  there  should  be  less  contraction,  just  as  we  had  found 
for  RBL  1  and  RBL  2.  Moreover,  the  performance  predictions 
are  in  even  poorer  agreement  with  the  test  data;  see  Table 
III. 


The  corresponding  fourth  approximation  to  the  x-varia¬ 
tions  of  the  induced  velocity  components  was  computed.  At 
the  0.4  and  0.7  stations,  the  computed  x-variations  of  u 
agree  well  with  the  assumed  variations.  At  0.9,  the  com¬ 
puted  variation  increases  more  rapidly  with  x  than  assumed 
but  has  a  maximum  and  subsequent  decrease.  At  0.975,  the 
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computed  x-variation  of  u  is  initially  less  rapid  but  con¬ 
tinues  to  increase  without  the  appearance  of  a  maximum  like 
the  one  that  occurred  for  the  original  assumed  variation. 

This  concluded  our  computations  for  this  case.  The 
results  and  trends  resemble  RBL  1  and  RBL  2  very  closely, 
in  that  an  apparent  tendency  toward  the  force- free  state  in 
the  x-variations  of  u  leads  to  a  less  favorable  comparison 
of  the  predicted  performance  with  test  data  as  well  as  a 
greater  discropency  between  the  assumed  and  computed  radial 
inflow. 

The  technique  that  we  developed  for  achieving  conver¬ 
gence  of  the  tangential  inflow  is  worth  mentioning  at  this 
point.  It  wa3  used  for  all  subsequent  cases  and  is  based 
upon  the  results  for  RBL  1  and  RBL  2,  that  the  computed 
tangential  inflow  differs  bv  only  a  small  amount  from 
-Nr/4*r.  At  each  step  of  the  iteration  we  determined  for 
each  radial  station  the  ratio  of  the  computed  value  to 
-Nr/4wr.  These  ratios  were  then  multiplied  by  the  next 
approximation  to  -Nr/4wr  to  determine  the  corresponding 
approximation  to  the  tangential  inflow.  This  procedure  led 
rapidly  to  convergence. 

Finally,  we  should  remark  on  the  correlation  of  the 
regular  part  of  the  axial  inflow  that  had  been  found  for 
the  RBL.  Here  it  did  not  prove  useful  because  there  were 
large  discrepancies  among  all  the  approximations,  especially 
between  the  initial  and  final  one.  We  attempted  to  general¬ 
ize  the  correlation  further  to  account  for  the  changes  in 
r  from  step  to  step  as  follows.  A  rough  approximation  to  an 
average  strength  of  the  derivative  of  r  can  be  found  by 
dividing  its  maximum  value  by  the  radial  distance  between 
the  maximum  and  the  tip.  This  distance  does  not  vary  much 
during  the  iteration.  Therefore,  we  tried  to  improve  the 
correlation  for  each  application  of  the  Biot-Savart  operator 
in  Eq.  (2)  by  dividing  the  regular  part  of  the  output  axial 
inflow  by  the  maximum  of  the  input  r  as  well  as  by  multiply¬ 
ing  it  by  the  input  axial  inflow  at  the  0.975  radial  station 
Unfortunately,  the  results  indicate  almost  equally  large 
discrepancies  among  the  iterations  here,  too.  We  did  use 
both  correlations  as  a  guide  later  in  the  early  stages  of 
iteration  when  aid  was  needed  most,  but  without  much  success 


3(109652)  Propeller,  16.0”  Pitch  Setting 

The  next  case  considered  was  the  16.0°  pitch  setting, 
as  measured  at  the  0.7  radial  station,  of  the  three-bladed , 
full-scale  Curtiss-Wright  3(109652)  propeller.  The  blade 
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design  characteristics  are  given  in  Appendix  II.  Computa¬ 
tions  on  the  basis  of  the  initial  wake  hypothesis  had  been 
made  for  this  case  previously;  see  Refs.  10  and  11.  The 
airfoil  sectional  data  appropriate  to  a  tip  speed  of  785 
feet  per  second  and  the  zeroth  approximations  were  deter¬ 
mined  as  usual  from  Refs.  6  and  12  and  the  generalized 
actuator  disk  results.  The  same  set  of  x-variations  of  the 
axial  velocity  u  was  selected  that  we  had  used  initially 
for  RBL  1,  RBL  2,  and  the  17.0°  sotting  of  the  120  AF 
propeller. 

Five  successive  approximations  to  r  and  the  inflow  were 
computed  to  achieve  convergence  of  the  axial  inflow  to  with¬ 
in  t2.lt,  except,  as  before,  at  the  0.2  radial  station.  To 
investigate  this  difficulty  inboard,  we  increased  the  hub 
radius  from  0.125  to  0.167  and  computed  two  further  approxi¬ 
mations,  which  did  not  improve  the  situation.  The  seventh 
approximations  to  r  and  the  axial  and  tangential  inflow  are 
presented  in  Table  IV  with  the  zeroth  approximations  and 
the  initial  wake  hypothesis  results  for  comparison.  The 
assumed  and  computed  radial  inflows  again  indicate  that 
there  should  be  less  contraction  than  we  have  specified. 

The  performance  is  plotted  in  Figure  14  with  the  initial 
wake  hypothesis  and  Curtiss-Wright  test  data  (Ref.  10)  for 
comparison.  It  should  be  noted  that  the  hub  radius  change 
led  to  performance  changes  of  less  than  *1.01. 

The  seventh  approximation  to  the  x-variations  of  the 
induced  velocity  components  along  the  trailing  vortex  sheets 
was  computed.  The  computed  x-variations  of  u  were  qualita¬ 
tively  similar  to  those  for  the  17.0°  setting  of  the  120  AF 
propeller  for  the  comparable  input  variations.  Maxima  were 
observed  at  both  the  0.9  and  the  0.975  radial  stations. 

In  view  of  the  unsuccessful  trend  that  we  had  found  in 
the  previous  cases  by  modifying  the  x-variations  of  u,  and 
because  the  indicated  changes  lead  in  the  same  direction  here, 
we  terminated  computations  for  this  case.  Instead,  we  pro¬ 
ceeded  to  cases  with  lower  pitch  settings. 


65  AF  Propeller,  12.6°  Pitch  Setting  -  Final  Computations 

We  returned  at  this  time  to  the  case  of  the  65  AF  pro¬ 
peller  at  a  pitch  setting  of  12.6°,  which  we  had  considered 
originially;  see  Chapter  4.  We  resumed  the  computation 
where  we  had  left  off,  using  the  x-variations  of  the  axial 
velocity  u  that  we  had  found  to  be  most  promising  at  that 
earlier  time. 
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Figure  14.  Comparison  of  Theoretical  and  Experimental  Performance;  3(109652 
Propeller,  nRp  ■  785  fps . 
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Ten  successive  approximations  to  r  and  the  inflow  were 
computed  to  achieve  axial  inflow  convergence  to  within  ^2.51 
everywhere  except  at  the  0.9  radial  station,  where  it  was 
8.11  different.  This  large  number  of  approximations  was 
necessary  because  the  axial  inflow  was  very  low  outboard 
(less  than  one-half  of  previous  cases),  and  the  results 
were  very  sensitive  to  small  changes.  The  tenth  approxi¬ 
mation  to  r  and  the  axial  and  tangential  inflow  are  presented 
in  Table  V  along  with  the  zeroth  approximations  and,  for  a 
further  comparison,  the  Knight  and  Hefner  values  according 
to  Ref.  19.  The  assumed  and  computed  radial  inflows  do  not 
compare  favorably;  see  Figure  15.  Contrasted  with  the 
previous  cases,  though,  the  indications  here  are  that  more 
contraction  is  needed  outboard.  The  predicted  performance 
is  again  much  too  optimistic  compared  with  the  data  acquired 
at  this  pitch  setting  in  our  cast  program.  In  particular, 
the  p indicted  and  measured  thrust  coefficients  are  0.1156 
and  0.1200,  respectively;  power  coefficients,  0.0377  and 
0.0520,  respectively;  and  figures  of  merit,  83%  and  631, 
respectively. 

The  tenth  approximation  to  the  x-variations  of  the 
induced  velocity  components  along  the  trailing  vortex  sheets 
was  computed.  At  the  0.4  and  0.7  radial  stations,  the  com¬ 
puted  x-variations  of  u  increase  slightly  less  rapidly,  as 
has  been  the  usual  case.  At  the  0.9  station,  however,  the 
computed  variation  is  in  agreement  initially  but  then  falls 
below  by  a  small  amount,  while  at  0.975  the  computed  varia¬ 
tion  increases  more  rapidly  initially  but  as  a  maximum  that 
is  about  301  lower  than  the  assumed  asymptotic  value.  We 
did  not  modify  the  x-variations  of  u  and  compute  further, 
but  extrapolations  based  on  the  previous  cases  show  that 
there  would  be  a  distinct  improvement  in  the  agreement  with 
performance  test  data  if  this  were  done. 

At  this  stage  of  our  investigation,  we  reexamined 
several  choices  of  basic  input  parameters  that  had  been 
selected  much  earlier.  This  was  appropriate  because  by  now 
we  had  achieved  convergent  inflow  iterations  in  a  number  of 
cases . 

First,  we  reduced  the  distance  between  the  basic 
elements  of  the  trailing  vortex  sheets  from  tt/6  to  ir/12. 

The  recomputed  axial  inflow  in  the  tenth  approximation  dif¬ 
fered  by  less  than  -2.51  everywhere  along  the  blade  and  was 
even  closer  near  the  tip.  This  is  within  what  we  considered 
to  be  convergence  for  this  case.  The  recomputed  performance 
parameters  differed  by  less  than  il.OI.  We  concluded  that 
our  basic  spacing  was  adequate. 
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Radial  Inflow 
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BLADE  CIRCULATION  AND  AXIAL  AND  TANGENTIAL  INFLOW 
65  AF  Propeller,  B0.7R„  »  12.6* 


O**1  -  Zeroth  Approxiaation 

K  5  H  -  Knight  and  Hefner  Approxiaation  (Ref.  19) 
10tn  .  Tenth  Approxiaation 


Next,  we  increased  the  number  of  Glauert  coefficients 
in  the  least- squares  fit  of  the  blade  circulation  distribu¬ 
tion  r  from  10  to  IS,  while  retaining  the  same  20  input 
values.  This  reduced  the  maximum  difference  between  input 
values  and  the  fit  from  -0.51  to  -0.21.  However,  there  was 
a  noticeable  introduction  of  "wiggles"  into  the  derivative 
of  r  and  thus  in  the  singular  part  of  the  axial  inflow. 

These  "wiggles"  arise  from  the  amplification  effects  of  the 
higher  order  Glauert  coefficients  that  we  discussed  in 
Chapter  2.  The  resultant  changes  in  the  tenth  approximation 
to  the  axial  inflow  are  less  than  -2.21  inboard  of  the  0.85 
radial  station  but  are  -5.6%,  +12.51,  and  -12.01  at  the  0.9, 
0.95,  and  0.975  radial  stations,  respectively.  The  perfor¬ 
mance  parameters  changed  by  less  than  0.31.  We  feel,  then, 
that  the  increases  at  some  radii  are  balanced  by  decreases 
at  others,  and  therefore  the  "wiggles"  are  indeed  extraneous. 
We  concluded  that  10  coefficients  had  been  a  wise  choice  but 
that  further  research  on  the  best  method  of  fitting  r  is 
desirable. 

Finally,  we  permitted  the  trailing  vortex  sheet  elements 
to  deform  immediately  behind  the  lifting  line  instead  of 
after  two  undeformed  initial  steps.  Although  the  influence 
functions  are  nonintegrable  in  a  rigorous  sense,  as  we  have 
pointed  out  in  Chapter  2,  the  numerical  integration  proceeds 
smoothly  nevertheless.  The  differences  in  the  tenth  approxi¬ 
mation  to  the  axial  inflow  were  less  than  -2.91  except  at 
0.975,  where  it  was  +4.71.  These  are  larger  differences 
than  we  have  generally  accepted  for  convergence,  although 
the  performance  parameters  differ  by  less  than  tl.01.  We 
concluded  that  the  axial  inflow  is  sufficiently  sensitive 
to  the  amount  of  initial  distortion  of  the  trail;  »g  vortex 
elements  that  an  improved  criterion,  based  on  physical 
considerations ,  should  be  investigated"! 

This  concluded  computations  for  this  case.  With  the 
success  achieved  at  this  relatively  low  pitch  setting,  we 
turned  to  two  further  examples  in  the  hope  of  generalizing 
the  results  obtained  here. 


120  AF  Propeller,  10.0°  Pitch  Setting 

The  next  case  is  the  10.0°  pitch  setting,  as  measured 
at  the  0.7  radial  station,  of  the  120  AF  propeller  operating 
at  a  tip  speed  of  700  feet  per  second.  The  instantaneous 
velocity  components  in  the  propeller  wake  had  been  measured 
for  this  case  in  the  experimental  portion  of  this  investi¬ 
gation  in  the  same  manner  as  for  the  8.2°  pitch  setting  of 
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Che  65  AF  propeller.  The  azimuthal  variations  of  the  axial 
and  radial  velocity  components  in  the  wake  are  presented  in 
Appendix  III,  while  full  details  of  the  experimental  inves¬ 
tigation  are  contained  in  Ref.  14. 

The  zeroth  approximations  to  the  blade  circulation  r 
and  the  axial  and  tangential  inflows  were  determined  as 
usual.  The  x-variations  of  the  axial  velocity  u  inboard  of 
the  0.8  radial  station  were  chosen  from  a  consensus  of  all 
the  earlier  computations,  while  outboard  of  0.8  they  were 
chosen  to  be  identical  to  those  for  the  12.6°  pitch  setting 
of  the  6S  AF  propeller  because  it  was  the  most  closely 
similar  case  that  we  have  treated. 

A  first  approximation  to  the  inflow  was  computed.  The 
axial  inflow  was  considerably  higher  than  the  input  over 
midblade  but  was  strongly  negative  outboard  of  the  0.93 
radial  station.  These  results  are  almost  identical  to  those 
described  in  Chapter  4  for  the  8.2°  setting  of  the  65  AF 
propeller.  In  that  case  we  had  changed  the  x-variations  of 
u  outboard  without  achieving  local  improvement  in  the  axial 
inflow,  a  result  which  was  borne  out  later  when  we  showed 
that  changes  in  the  x-variations  of  u  outboard  do  not  change 
the  axial  inflow  significantly  outboard  (see  Figures  9,  11, 
and  13).  Accordingly,  we  retained  the  same  x-variations  of 
u  here,  but  we  chose  a  first  approximation  to  the  axial 
inflow  that  was  qualitatively  similar  to  our  earlier  con¬ 
verged  results.  Unfortunately,  the  second  approximation  was 
even  acre  strongly  negative  and  showed  even  less  change  of 
achieving  inflow  convergence. 

The  wake  data  were  again  limited  in  their  usefulness  in 
exactly  the  samo  way  that  they  were  for  the  65  AF  propeller 
at  8.2°;  see  Chapter  4.  The  extent  and  growth  downstream  of 
the  nonrepeatable  region  are  presented  in  Figure  8. 

Therefore,  this  case  as  well  as  the  earlier  one  could 
not  be  iterated  further  with  our  present  knowledge.  A  good 
deal  of  further  research,  both  experimental  and  theoretical, 
is  required  to  understand  the  nature  of  the  flow  phenomena 
in  the  vicinity  of  the  propeller  tip,  especially  in  cases 
»uch  as  these. 


3(109652)  Propeller,  10.0°  Pitch  Setting 

The  final  case  considered  is  the  10.0°  pitch  setting, 
as  measured  at  the  0.7  radial  station,  of  the  3(109652) 
propeller  operating  at  a  tip  speed  of  785  feet  per  second. 
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This  case  had  been  treated  in  our  earlier  studies  based  upon 
the  initial  wake  hypothesis,  and  a  single  approximation  hr.d 
also  been  calculated  with  the  refined  wake  hypothesis.  The 
zeroth  approximations  to  r  and  the  axial  and  tangential  in¬ 
flows  were  found  as  usual,  while  the  same  set  of  x-vari&tions 
of  u  was  selected  that  we  had  tried  at  the  10.0*  pitch  set¬ 
ting  of  the  120  AF  propeller. 

Eight  successive  approximations  to  r  and  the  inflow 
were  computed.  The  axial  inflow  converged  inboard  to  with¬ 
in  *2.01,  although  outboard  the  differences  were  as  large  as 
-6.1%  despite  the  absolute  differences  being  of  the  saae 
size.  The  axial  inflow  near  the  tip  is  the  lowest  of  any 
case  we  have  computed.  Again,  convergence  could  not  be 
achieved  at  the  0.2  radial  station.  The  eighth  approxi¬ 
mations  to  r  and  the  axial  and  tangential  inflows  are  pre¬ 
sented  in  Table  VI  along  with  the  zeroth  approximations  and 
the  initial  wake  hypothesis  results.  The  assumed  and  com¬ 
puted  radial  inflows  do  not  compare  well,  and  the  indications 
are  that  more  contraction  is  required  outboard,  as  was  shown 
in  Figure  15  for  the  65  AF  propeller  at  12.6*.  This  is  in 
marked  contrast  to  the  16. 0*  setting  of  this  same  3(109652) 
propeller.  The  predicted  performance  is  again  too  optimistic 
compared  to  test  data;  see  Figure  14. 

The  eighth  approximation  to  the  x-variations  of  the  in¬ 
duced  velocity  components  along  the  trailing  vortex  sheets 
was  computed.  At  the  0.4  and  0.7  radial  stations,  the  com¬ 
puted  x-variations  of  u  lie  slightly  below  the  assumed  values. 
At  the  0.9  station,  the  computed  variation  increases  more 
rapidly,  has  a  maximum  which  is  larger  than  the  assumed 
asymptotic  value,  and  then  decreases.  At  0.975,  the  computed 
variation  increases  initially  at  the  assumed  rate  but  has  a 
maximum  that  is  less  than  one-half  as  large  as  the  assumed 
asymptote.  Qualitatively,  the  indications  for  this  case  are 
the  same  as  for  the  65  AF  propeller  at  12.6*;  namely,  that 
lower  asymptotic  values  of  the  x-variations  of  u  are  needed 
outboard.  Extrapolations  based  on  such  modified  variations 
should  give  performance  predictions  that  are  in  better  agree¬ 
ment  with  the  test  data.  These  extrapolations  have  not  been 
verified . 


Observations  for  Specific  Propellers 

The  additional  complications  of  the  blade  circulation 
r  changing  from  one  approximation  to  the  next,  as  well  as 
the  geometrical  differences  in  blade  design  from  one  propel 
ler  to  another,  limit  our  observations  and  generalizations 


BLADE  CIRCULATION  AND  AXIAL  AND  TANGENTIAL  INFLOW 
3(109o.2)  Propeller,  80.7R-  -  10.0* 


Zeroth  Approximation 
Eighth  Approximation 
Initial  Wake  Hypothesis  (Ref.  10) 


to  a  greater  extent  than  for  the  representative  blade 
loadings.  Moreover,  we  have  already  made  observations  about 
our  attempt  to  generalize  the  correlation  of  the  regular 
part  in  the  section  for  the  120  AP,  17.0*  pitch  setting  case. 

The  first  area  of  observations  concerns  the  inflow 
iterations  for  a  given  set  of  x-variations  of  u.  Unfortu¬ 
nately.  the  iterative  properties  for  the  axial  inflow  have 
not  emerged  clearly.  The  sensitivity  that  we  discussed 
earlier  regarding  the  fit  of  r  contributes  to  the  difficulty. 
For  instance,  the  convergence  did  begin  outboard  generally 
and  proceeded  inboard.  However,  this  trend  is  not  as 
decisive  as  for  the  representative  blade  loadings  because 
r  changed  as  well,  and  so  the  changes  in  the  axial  inflow 
outboard  from  step  to  step  were  larger  than  the  corresponding 
changes  for  a  fixed  r.  Nevertheless,  the  changes  were 
definitely  smal’er  outboard  than  inboard.  Also,  it  was  much 
more  difficult  to  get  the  axial  inflow  outboard  even  to 
settle  down  enough  to  continue  the  iteration;  we  could  not, 
in  fact,  for  the  two  cases  for  which  we  obtained  measured 
wake  data.  We  are  unable  to  recommend  precise  values  for 
the  axial  iteration  factors,  but  we  found  that  values 
between  0.5  and  1.0  were  best  inboard  of  the  0.8S  radial 
station  and  that  values  between  0.25  and  0.5  were  best 
outboard  of  0.85.  A  successful  technique  for  iterating  the 
tangential  inflow  has  been  discussed  in  the  section  for  the 
120  AF  propeller  at  the  17.0°  pitch  setting. 

In  several  cases,  a  convergent  axial  inflow  could  not 
be  achieved  at  the  0.2  radial  station,  due  in  part  to  the 
blade  section's  being  beyond  the  angle  of  attack  of  normal 
stall.  In  view  of  tne  small  effect  that  the  inboard  region 
has  on  the  performance,  as  witnessed  by  the  small  performance 
change  with  the  hub  radius  modification  for  the  3(109652) 
propeller  at  the  16.0°  pitch  setting,  it  would  be  sensible 
just  to  fix  the  hub  radius  consistently  at  the  0.2  station 
in  all  cases. 

The  second  area  of  observations  concerns  the  iteration 
toward  a  satisfactory  force-free  approximation  by  means  of 
changes  in  the  x-variations  of  u.  The  single  case  that  we 
did  iterate  in  this  fashion  behaved  identically  to  RBL  1 
and  RBL  2;  namely,  the  x-variations  of  u  seemed  to  be 
approaching  the  force- free  state,  but  the  computed  radial 
inflow  indicated  r  en  more  strongly  that  less  contraction 
was  needed  than  had  been  assumed.  Furthermore,  the  agree¬ 
ment  between  performance  predictions  and  test  data  became 
worse.  This  case  had  a  thrust  coefficient  higher  than  both 
RBLs .  Unfortunately,  no  changes  were  made  in  the  x-variations 
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of  u  after  the  first  inflow  itoration  for  the  two  lower 
pitch-setting,  lower  thrust-coefficient  cases.  For  these, 
however,  the  computed  radial  inflow  indicated  that  more 
contraction  was  needed  than  had  been  assumed.  It  is  clear, 
then,  that  several  problems  must  be  resolved  before  a 
sufficiently  good  approximation  to  the  performance  can  be 
achieved. 
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CONCLUSIONS 


Based  on  the  results  obtained  in  this  study,  we  draw 
the  following  principal  conclusion. 

Instantaneous  wake  velocity  data  that  were  repeatable 
with  each  blade  passage  could  be  obtained  only  over  the 
middle  part  of  the  blade  for  the  65  AF  propeller  at  the  8.2* 
pitch  setting  and  for  the  120  AF  propeller  at  the  10.0°  set¬ 
ting.  Near  the  hub  and  especially  near  the  tip,  the  data 
were  random  and  could  not  be  resolved  meaningfully  into 
components.  Theoretical  computations  for  these  two  cases, 
which  are  low  pitch  settings,  could  not  be  iterated  success¬ 
fully.  Lack  of  data  near  the  tip  prevented  guidance  of  the 
computations,  and  lack  of  convergence  prevented  checking  with 
the  good  data. 

The  refined  wake  hypothesis  was  found  to  be  inadequate, 
in  general,  as  an  approximation  to  the  force-free  condition 
of  the  wake  on  the  basis  of  poor  agreement  of  performance 
predictions  with  test  data  and  an  excessive  amount  of  con¬ 
traction.  Accordingly,  a  new  wake  representation  was  for¬ 
mulated  which  fixed  the  overall  contraction  pattern  and 
permitted  concentration  on  the  effective  "pitch"  of  the 
trailing  vortex  sheets. 

A  series  of  computations  for  thrust  coefficients  above 
0.15  (namely,  two  representative  blade  loadings  and  two 
propellers,  each  at  one  setting)  showed  that  the  assumed 
and  computed  "pitch"  came  into  closer  agreement  (i.e.,  more 
nearly  force  free)  through  changes  in  the  variations  of  the 
axial  velocity  with  distance  downstream  along  the  trailing 
vortex  sheets.  Simultaneously,  though,  the  fixed  contraction 
became  progressively  too  great  near  the  tip,  and  the  perfor¬ 
mance  predictions  became  too  optimistic  compared  to  test 
data . 


A  series  of  computations  for  thrust  coefficients  below 
0.13  (namely,  two  propellers,  each  at  one  setting)  was  made 
without  changing  the  variations  of  the  axial  velocity.  The 
fixed  contraction  for  these  cases  was  too  small  near  the 
tip,  but  the  performance  predictions  were  still  too  optimis¬ 
tic. 
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Practical  use  as  a  design  tool  of  solutions  for 
specific  blade  circulation  distributions,  such  as  the 
representative  blade  loadings,  is  premature  without  better 
theoretical  and  experimental  correlation.  However,  the 
representative  loadings  offer  great  computational  simplifi¬ 
cations  that  can  be  utilized  in  future  theoretical  develop¬ 
ments. 
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APPENDIX  I 


OPERATING  INSTRUCTIONS  FOR  COMPUTER  PROGRAMS 


General 


Brief  operating  instructions  are  presented  to  outline 
the  use  of  the  computer  programs  for  the  inflow  velocity 
components  (INFLOWPT),  for  the  induced  velocity  components 
along  the  trailing  vortex  sheets  (SHEETPT) ,  for  the  propel¬ 
ler  blade  circulation  and  performance  (PROPERFM) ,  and  for 
the  vortex  element  trajectories  and  associated  influence 
functions  (MRCHINFO) .  These  instructions  consist  of  descrip¬ 
tions  of  the  input  data  required,  the  punch  format,  and  the 
print-out.  All  angles  are  in  degrees,  and  the  remaining 
data,  unless  noted  otherwise,  must  be  nondimensionalized 
with  respect  to  the  characteristic  time  1/n  and  the  charac¬ 
teristic  length  Rp. 


Program  INFLOWPT 

It  is  assumed  that  the  following  basic  data  are  known; 
namely,  distributions  of  the  blade  circulation  and  inflow 
as  well  as  the  variations  with  axial  distance  downstream  of 
the  induced  velocity  components  along  the  trailing  vortex 
sheets. 

Specify  the  number  of  propeller  blades  (NB) ,  the  pro¬ 
peller  hub  radius  (RH) ,  the  total  number  of  points  at  which 
the  inflow  is  desired  (NI),  and  the  radial  coordinates  (RI) 
of  these  inflow  points.  NI  cannot  exceed  40. 

Specify  the  number  of  Glauert  coefficients  (NC)  and 
the  input  values  of  the  blade  circulation  (G)  at  each  of  the 
total  number  (NG)  of  Glauert  variable  locations  (PG) . 

NC  and  NG  cannot  exceed  25  and  100,  respectively. 

Specify  the  normalized  variations  along  the  trajectories 
of  the  axial  (UW)  and  tangential  (WW)  velocity  components 
at  a  total  number  (JW)  of  radial  stations  (RW)  for  each  of 
a  total  number  (IW)  of  axial  stations  (XW).  Specify  the 
ratio  of  radial -to-axial  velocity  (VW)  along  the  trajectories 
for  the  same  stations.  Specify  the  magnitudes  of  the  free- 
.itream  velocity  (UFS)  and  the  axial  (UFAC)  and  tangential 
(WFAC)  inflow  components  for  the  same  JW  radial  stations. 

IW  and  JW  cannot  exceed  30  each.  In  the  present  form  of 
the  program,  VFAC  is  not  used,  but  zeroes  must  be  read  in. 
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Specify  the  number  of  steps  (NSTEPSTR)  and  the  associ¬ 
ated  step  size  of  the  time  increment  (STBPSTRT)  for  the  ini¬ 
tial  undeformed  march  from  the  blade.  For  the  remainder  of 
the  march,  specify  the  number  of  steps  (NSTBP)  at  each  of 
the  total  number  (NINCSIZE)  of  marching  step  time  Increment 
sizes  (STEP).  NINCSIZE  cannot  exceed  10.  Specify  the 
axial  location  (XCUTOPF)  at  which  the  march  is  terminated; 
this  location  must  coincide  with  XW(IW).  Specify  the  abso¬ 
lute  value  of  accuracy  (ERR)  required  in  computing  the 
Heuman  lambda  function  (0.0001  was  used).  Specify  the 
parameter  (NE)  that  determines  the  fraction  of  *,  i.e., 
w/NE,  between  the  basic  trailing  vortex  elements. 

Specify  the  absolute  value  of  accuracy  (ACC)  required 
in  the  Gaussian  integration  scheme  and  the  maximum  number 
of  subdivisions  (NLIMIT)  permitted  before  termination. 

Specify,  finally,  the  month  (IM),  day  (ID),  aind  year 
(IY)  of  the  run,  the  run  number  (NUMRUN) ,  and  a  title  for 
the  run  (TITLE).  TITLE  cannot  exceed  80  alphanumeric 
characters . 

The  punch  format  for  program  INFLOWPT  is  given  in 
Table  VII  and  includes  the  input  quantities  and  the  field 
specification. 

The  print-out  of  INFLOWPT  consists  of  the  following: 
a  title  page;  the  axial,  radial,  and  tangential  velocity 
components  of  the  stored  input;  identification  of  the 
automatically  selected  elements  of  the  trailing  vortex 
sheets;  the  Glauert  coefficients  and  the  fit  of  the  blade 
circulation  and  its  derivative  with  respect  to  the  Glauert 
variable;  and  the  strengths  of  the  Cauchy  and  logarithmic 
singularities  as  defined  in  Eqs.  (24)  through  (26).  For 
each  trailing  element  the  following  are  printed  out  at  all 
inflow  points:  the  contributions  to  the  total  influence 
functions  from  the  deformed  portion  of  the  trajectory;  the 
total  influence  functions;  and  the  regular  parts  of  the 
total  influence  functions.  For  each  Inflow  point  there  is 
printed  out  the  regular  part  of  the  total  influence  functions 
from  all  trailing  elements  contributing  to  that  point.  The 
final  results  consist  of  the  following  print-out:  the  con¬ 
tribution  to  the  inflow  components  from  the  regular  part  of 
the  integrand;  the  contribution  from  the  Cauchy  singular 
part;  the  contribution  from  the  logarithmically  singular 
part;  the  total  singular  contribution;  and  the  total  inflow 
components . 
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TABLE  VII 


INPUT  PUNCH  FORMAT  •  PROGRAM  INFLOWPT 


Inyut  Quintity 
IM,  ID,  IY 
NUMRUN 

TITLE(l)  -  TITLE(IO) 

NB,  NE,  NI,  NG,  NC,  NINCSIZE 

NLIMIT,  IN,  JW 

RH,  XCUTOPF 

ACC,  ERR 

XW(1)  -  XW(IW) 

RW(1)  -  RW (JW) 

UW(1,1)  -  UW(IW,1) 

UWCI.JW)’-  UWJW.JW^ 

VW(1,1)  -  VW(IW,1) 

vw(i,jw)’-  vw(iw,.jw2_ 

WW(1,1)  -  WW(IW,1) 

« 

mam 

WW(I,JW)’-  WW( IN , JW) 


Field 
Spec . 

16 1  S 

161  S 

1  0A8 

1615 

16  I  S 

10F8.5 

10F8.5 

10F8.5 

10F8.5 

10F8.5 


10F8.5 


10F8.S 
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TABLE  VH-contd. 


Program  SHEETPT 

Tho  same  basic  data  are  required  as  for  INFLOWPT. 

Specify  the  number  of  propeller  blados  (NB) ,  the  pro- 
pellor  hub  radius  (RH)  ,  the  numbor  of  axial  positions  down¬ 
stream  (NXF)  at  which  the  sheet  velocities  are  to  be  com¬ 
puted  and  their  locations  (XF)  ,  the  number  of  trajectories 
(NRB)  along  which  tho  shoot  velocities  aro  to  bo  computed 
and  thoir  radial  coordinates  at  the  blade  (RB) ,  and  the 
number  of  added  auxiliary  trajectories  (NRA)  adjacent  to  th? 
sheet  velocity  trajectories  and  their  radial  coordinates  at 
the  blade  ( RA ) .  NXF,  NRB,  and  NRA  cannot  exceed  10,  20, 
and  40,  respectively. 

Specify  the  number  of  Glauert  coefficients  (NC)  and 
the  input  values  of  the  blado  circulation  (G)  at  each  of 
the  total  number  (NG)  of  Glauert  variable  locations  (PG). 
Neither  NC  nor  NG  can  exceed  100. 

Specify  the  normalized  variations  along  the  trajectories 
of  the  axial  (UW)  and  tangential  (WW)  velocity  components 
at  a  total  number  (JW)  of  radial  stations  (RW)  for  each  of 
a  total  number  (IW)  of  axial  stations  (XW).  Specify  the 
ratio  of  radial-to-axial  velocity  (VW)  along  the  trajectories 
for  the  same  stations.  Specify  the  magnitudes  of  the  free- 
stream  velocity  (UFS),  and  the  axial  (UFAC)  and  tangential 
(WFAC)  inflow  components  for  the  same  JW  radial  stations. 

IW  and  JW  cannot  exceed  30  each.  In  the  present  form  of  the 
program,  VFAC  is  not  used,  but  zeroes  must  be  read  in. 

Specify  the  number  of  steps  (NSTEPSTR)  and  the  associ¬ 
ated  step  size  of  the  time  increment  (STEPSTRT)  for  the  ini¬ 
tial  undeformed  march  from  the  blade.  For  the  remainder  of 
the  march  up  to  XF(1)  and  for  the  march  after  each  XF, 
specify  the  number  of  steps  (NSTEP)  at  each  of  the  total 
number  (NINCSIZE)  of  marching  step  time  increment  sizes 
(STEP).  NINCSIZE  cannot  exceed  10.  Specify  the  marching 
step  time  increment  size  (STEPT)  that  is  used  within  an 
axial  interval  (DELTAX)  on  either  side  of  each  XF.  Specify 
the  axial  location  (XCUTOFF)  at  which  the  march  is  terminated; 
this  location  must  coincide  with  XW(IW).  Specify  the  abso¬ 
lute  value  of  accuracy  (ERR)  required  in  computing  the  Heuman 
lambda  function  (0.0001  was  used).  Specify  the  parameter 
(NE)  that  determines  the  fraction  of  w ,  i.e.,  w/NE,  between 
the  basic  trailing  vortex  elements. 
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Specify  the  absolute  value  of  accuracy  (ACC)  required 
in  the  Gaussian  integration  scheme  and  the  maximum  number 
of  subdivisions  (Nl.IMIT)  permitted  before  termination. 

Use  the  following  values:  for  (RLIMT),  0.987;  for 
(RVDNA) ,  0.02S;  for  (RVDNB )  ,  1.950;  for  (RVUPA)  ,  2.050; 
for  (RVUPB) ,  2.040;  and  for  (RVUPC),  2.026. 

Specify,  finally,  the  month  (IM),  duy  (ID,  and  year 
(1Y)  of  the  run,  the  run  number  (NUMRUN) ,  and  u  title  for 
the  run  (TITLE).  TITLE  cannot  exceed  80  alphanumeric 
characters . 

The  punch  format  for  program  SHEETPT  is  given  in  Table 
VIII  and  includes  the  input  quantities  and  the  field  speci- 
f icat ion . 

The  print-out  of  SHEETPT  consists  of  the  following:  a 
title  page;  the  axial,  radial,  and  tangential  velocity  com¬ 
ponents  of  the  stored  input;  identification  of  the  auto¬ 
matically  selected  elements  of  the  trailing  vortex  sheets; 
the  Glaucrt  coefficients  and  the  fit  of  the  blade  circulation 
and  its  derivative  with  respect  to  tho  Glauert  variable; 
identification  of  the  coordinates  of  the  auxiliary  points 
along  the  trailing  vortex  sheets  that  arc  used  in  determining 
y  and  the  derivative  of  S;  identification  of  the  coordi¬ 
nates  of  the  sheet  points;  and  the  strengths  of  the  Cauchy 
and  logarithmic  singularities  as  defined  in  Table  II  along 
with  V  and  the  derivative  of  S.  For  each  trailing  element 
the  following  are  printed  out  at  all  sheet  points:  the 
contributions  to  the  total  influence  functions  from  the 
deformed  portion  of  the  trajectory;  the  total  influence 
functions;  and  the  regular  parts  of  the  total  influence 
functions.  For  each  sheet  point  there  is  printed  out  the 
regular  part  of  the  total  influence  functions  from  all 
trailing  elements  cintributing  to  that  point.  The  final 
results  consist  of  the  following  print-out:  the  contribution 
to  the  sheet  velocity  components  from  the  regular  part  of 
the  trailing  vortex  integrand;  the  contribution  from  the 
Cauchy  singular  part;  the  contribution  from  the  logarithmic¬ 
ally  singular  part;  the  total  singular  contribution;  the 
total  trailing  vortex  sheet  contribution  (i.c.,  the  sum  of 
the  regular  and  singular  contributions);  the  contribution 
from  the  bound  blade  vortices;  and  the  total  induced  velocity 
components  at  the  sheet  points. 


TABLE  VIII 


INPUT  PUNCH  FORMAT  -  PROGRAM  SHEETPT 


Input  Quantity 

Field 
Spec . 

IM,  ID,  IY 

1615 

NUMRUN 

1615 

TITLE  (1)  -  TITLE  (10) 

1 0A8 

NB,  NE,  NRB,  NRA,  NXF ,  NG ,  NC 

1615 

NLIMIT,  IW,  JW 

1615 

RH,  XCUTOFF ,  DELTAX, 

STEPT 

10F8.5 

ACC,  ERR 

10F8.5 

XW(1)  *  XW(IW) 

10F8.5 

RW(1)  -  RW (JW) 

10F8.5 

UW(1,1)  •  UW(IW,  1)" 

• 

• 

10F8.5 

UW(1,JW)*-  UW(IW,JW) 

VW(1,I)  -  VW(IW,1) 

• 

• 

10F8.5 

VW(1,JW)’-  VW(IW,JW) 

WW(1,1)  -  WW(IW,1) 

• 

• 

• 

! 

10F8.5 

WW(1,JW)  -  WW(IW,JW) 


TABLE  VIII -contd . 


UPS ( 1 )  -  UFS(JW) 

UFAC(l)  -  UFAC(JW) 

VFAC(l)  -  VFAC(JW) 

WFAC(l)  -  WFAC(JW) 

RB ( 1 )  -  RB(NRB) 

RA(1 )  -  RA(NRA) 

XF  ( 1 )  -  XF(NXF) 

PG(1)  -  PG(NG) 

G C 1 )  -  G(NG) 

NSTEPSTR,  STEPSTRT 
NINCSIZE(l)  “ 

NSTEP (1,1)  -  NSTEP(1,NINCSIZE(1)) 
STEP (1,1)  -  STEP(1,NINCSIZE(1)) 


10P8.S 
10F8.S 
10F8.5 
10F8.5 
10F8.5 
10F8.5 
10F8.5 
10F8.5 
10F8.S 
I  5 , F8 . 5 
1615 
1615 
10F8.5 


NINCSI ZE  (NXF* 1 ) 

NSTEP  (NXF*  1  , 1 )  -  NSTEP  (NXF«- 1  ,NINCSIZE(NXF*1) 
STEP (NXF* 1 , 1 )  -  STEP(NXF*1  ,NINCSIZE (NXF*1 ) ) 
RLIMT ,  RVDNA,  RVDNB ,  RVUPA,  RVUPB,  RVUPC 


1615 

1615 

10F8.5 


10F8.5 


Program  PROPERFM 

It  is  assumed  that  the  following  basic  data  are  known; 
namoly,  the  geometrical  characteristics  of  the  propeller, 
the  sectional  lift  and  drag  data,  and  the  inflow. 

Specify  the  number  of  propeller  blades  (NB) ,  the  hub 
radius  (RIIB) ,  the  tip  radius  in  feet  (RP) ,  and  the  tip 
speed  in  feet  per  second  (PTS). 

Specify  at  each  of  the  total  number  (NRB)  of  radial 
stations  (RB)  the  blade  chord  in  feet  (B) ;  the  Mach  number 
(MN) ;  the  Reynolds  number  times  ten  to  the  minus  sixth 
power  (REN);  a  quantity  (ICDCURVE)  that  is  1  if  the  drag 
coefficient  (YCD)  at  that  RB  is  given  as  a  function  of  angle 
of  attack  (XAL)  and  is  2  if  YCD  is  given  as  a  function  of 
lift  coefficient  (XCL)  ;  a  quantity  (CDZ)  which  is  the  drag 
coefficient  for  a  20°  angle  of  attack  at  that  RB  if  ICDCURVE 
■  1  or  is  the  drag  coefficient  for  an  8°  angle  of  attack  if 
ICDCURVB  ■  2.  NRB  cannot  exceed  15. 

Specify  the  total  number  of  points  (NXAL)  at  which  XAL 
is  given  and  the  total  number  of  points  (NXCL)  at  which  XCL 
is  given.  Specify  for  each  RB  the  lift  coefficient  (YCL) 
at  each  value  of  XAL,  and  specify  YCD  at  each  value  of  XAL 
if  ICDCURVE  -  1  or  at  each  value  of  XCL  if  ICDCURVE  -  2. 
Specify  a  factor  (CDFAC)  to  multiply  the  drag  coefficient, 
if  desired;  usually  1.0  was  taken.  Neither  NXAL  nor  NXCL 
can  exceed  60. 

Specify  for  each  RB  the  free-stream  velocity  (CUB), 
the  axial  (SUB)  and  tangential  (SWB)  inflow,  and  the  blade 
pitch  angle  (BETA) . 

Specify,  finally,  the  month  (IM),  day  (ID),  and  year 
(IY)  of  the  run,  the  run  number  (NUMRUN)  and  case  (NUMCASE)  , 
and  a  title  for  the  case  (PROPID).  PROPID  cannot  exceed  80 
alphanumeric  characters. 

The  punch  format  for  program  PROPERFM  is  given  in  Table 
IX  and  includes  the  input  quantities  and  the  field  specifi¬ 
cation.  Note  that  the  YCD  input  is  given  for  ICDCURVE  -  1; 
if  ICDCURVE  ■  2,  the  upper  limit  would  be  NXCL  rather  than 
NXAL. 


The  print-out  of  PROPERFM  consists  of  a  title  page 
first.  Tnen,  for  each  radial  station,  all  the  input 
quantities  described  above  are  listed,  including  the  lift 
and  drag  data  as  well  as  the  following  computed  quantities: 
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TABLE  IX 


INPUT  PUNCH  FORMAT  -  PROGRAM  PROPERFM 


Field 


Input  Quantity 

Spec. 

IM,  ID,  IY 

1615 

NUMRUN,  NUMCASE 

1615 

NB,  NRB,  NXAL,  NXCL 

1615 

RP,  RHB ,  PTS 

10F8.5 

RB Cl)  -  RB(NRB) 

10F8.5 

B (1)  -  B(NRB) 

10F8.5 

CDZ(l)  -  CDZ(NRB) 

10F8.5 

MN(1)  -  MN(NRB) 

10F8.5 

REN(l)  -  REN (NRB) 

10F8.5 

XAL(l)  -  XAL(NXAL) 

10F8.5 

XCL(l)  -  XCL(NXCL) 

10F8.5 

ICDCURVE (1)  -  ICDCURVE(NRB) 

1615 

YCL  (1 , 1)  -  YCL (1 ,NXAL) 

10F8.5 

j 

YCD(l.l)  -  YCD (1 ,NXAL) 

• 

10F8.5 

I 

• 

YCL (NRB, lj  -  YCL (NRB, NXAL  fl 

10F8.5 

YCD (NRB, 1)  -  YCD (NRB, NXAL  )J 

mm 

I0F8.5 

PROPID(l)  -  PROPID (10) 

10A8 

- .  ^  a.  ' 


TABLE  IX-contd. 


CUB(l)  -  CUB(NRB) 
SUB(l)  -  SUB(NRB) 
SWB(l)  -  SWB(NRB) 
BETA (1)  -  BETA(NRB) 
CDFAC 


10F8.5 

10F8.5 

10F8.5 

10F8.5 


10F8.5 
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the  Glauert  variable  of  the  radial  station,  the  nondimen- 
sionalized  blade  chord,  the  total  local  velocity  seen  by 
the  station,  the  angle  between  the  rotational  velocity  and 
total  velocity,  the  local  blade  angle  of  attack,  the  local 
lift  and  drag  coefficients,  the  local  blade  circulation  r, 
the  nondimensionalized  value  of  -Nr/4irr  from  the  local  r, 
and  the  radial  derivatives  of  the  thrust  and  power  coeffi¬ 
cients.  On  a  first  summary  page,  all  the  above  computed 
quantities  are  repeated,  and  the  integrated  thrust  and  power 
coefficients  and  figure  of  merit  are  presented.  On  a  second 
summary  page,  the  data  on  the  first  summary  page  are  repeated 
but  with  the  drag  coefficients  all  set  identically  equal  to 
zero.  Finally,  the  integrated  induced  thrust  and  power 
coefficients  and  figure  of  merit  are  presented. 


Program  MRCHINFO 

This  program  was  abstracted  from  INFLOWPT  and  assumes 
knowledge  of  most  of  the  same  basic  data.  There  is  some 
flexibility  in  the  use  of  MRCHINFO,  but  we  have  selected 
the  mode  that  gives  the  total  influence  functions  at  a 
number  of  inflow  points  as  well  as  the  trajectory. 

Specify  (IM)  ,  (ID),  (IY),  (NUMRUN) ,  (TITLE),  (IW),  (JW) , 
(XW),  (RW),  (UW),  (VW),  (WW),  (UFS) ,  (UFAC) ,  (VFAC) ,  (WFAC)  , 
(RH),  (NB),  (NI),  (RI),  (NSTEPSTR),  (NINCSIEE)  ,  (STEPSTRT)  , 
(XCUTOFF) ,  (ERR),  (NSTEP)  and  (STEP)  exactly  the  same  as 
for  program  INFLOWPT.  Specify  further  a  case  number  (ICASE) . 
Use  the  following  values:  for  (INFO),  2;  for  (IDATA) ,  1. 
Specify  the  coordinate  of  the  desired  trailing  vortex  element 
(VORTEXPT) .  This  coordinate  may  be  expressed  in  blade  radii 
or  in  the  Glauert  variable  in  degrees,  provided  the  latter  is 
greater  than  1°.  Specify  a  factor  (CONTRFAC)  that  multiplies 
VW  uniformly  to  provide  for  a  different  amount  of  contraction 
than  that  fixed  by  the  basic  VW  field. 

The  punch  format  for  program  MRCHINFO  is  given  in 
Table  X  and  includes  the  input  quantities  and  the  field 
specification. 

The  print-out  of  MRCHINFO  consists  of  the  following: 
a  title  page;  the  axial,  radial,  and  tangential  velocity 
components  of  the  stored  input;  and  identification  of  the 
number  of  steps  at  each  of  the  marching  step  time  increment 
sizes.  For  the  trailing  element  the  following  are  printed 
out:  at  each  step  of  the  march,  the  coordinates  and 
velocity  components  together  with  the  cumulative  contribu¬ 
tions  to  the  axial,  radial,  and  tangential  influence 
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TABLE  X 


INPUT  PUNCH  FORMAT- PROGRAM  MRCHINFO 


Field 

Input  Quantity 

Spec . 

IM,  ID,  IY 

1615 

INFO,  I DATA 

1615 

NUMRUN,  I CASE 

1615 

TITLE(l)  -  TITLE (10) 

10A8 

IW,  JW 

1615 

XW(1)  -  XW(IW) 

10F8.5 

RW(1)  -  RW (JW) 

10F8.5 

UW(1,1)  -  UW(IW,1) 

• 

• 

—  10F8.5 

UW(1,JW)*-  UW(IW,JW)_ 

VW,(1,1)  -  VW(IW,1)  " 

• 

• 

VW(1,JW)  -  VW(IW,JW)_ 

—  10F8.5 

WW(1,1)  -  WW(IW,1) 

• 

• 

-  10F8.5 

WW(1,JW)  -  WW(IW,JW) 

UFS(l)  -  UFS(JW) 

10F8.5 

UFAC(l)  -  UFAC(JW) 

10F8.5 

j 

| 
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TABLE  X-contd. 


VFAC(l)  -  VFAC(JW)  10F8.5 
WFAC(l)  -  WFAC(JW)  10F8.5 
CONTRFAC  10F8.S 
VORTEXPT,  RH  10F8.5 
NB.NI  1615 
R I  ( 1 )  -  RI(NI)  10F8.5 
NSTEPSTR ,  NINCSIZE,  STEPSTRT , XCUTOFF  ,  ERR  2I5.3F8.5 
NS1EP(1)  -  NSTEP  (NINCSIZE)  1615 
STEP ( 1 )  -  STEP (NINCSIZE)  10F8.5 
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functions  at  all  the  Inflow  points;  the  contributions  to 
the  total  influence  functions  from  the  deformed  portion  of 
the  trajectory;  and  the  total  influence  functions. 


APPENDIX  II 


BLADE  CHARACTERISTICS 

Three  propellers  are  considered  in  the  text  for  both 
theoretical  calculations  and  comparison  with  test  data; 
namely,  the  four-bladed  Canadair  65  Activity  Factor  and  120 
Activity  Factor  propellers  and  the  three-bladed  Curtiss- 
Wright  3(109652)  propeller.  The  digits  within  the  parenthe¬ 
ses  of  the  latter  refer  to  the  blade  design.  The  geometries 
and  characteristics  of  these  blade  designs  are  given  in 
Figures  16  through  18,  with  h/b  as  the  local  blade  thickness- 
to- chord  ratio  and 


AF 


100,000 

16 


b 

2RP 


(27) 


(28) 


where 

AF  is  the  activity  factor  of  the  blade. 


ICj^.  is  the  integrated  design  lift  coefficient  of  the 
1  blade,  and 

C,  is  the  blade  sectional  design  lift  coefficient 
^i  (Ref.  6,  pp.  2-5,  for  example). 
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gure  17.  Blad 
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Figure  18.  Blade  Characteristics,  109652  Design,  NACA  65-Series  Sec 
AF  =  115,  IC,  =  0.500. 


APPENDIX  III 


RESULTS  OF  INSTANTANEOUS  WAKE  VELOCITY  MEASUREMENTS 

In  the  experimental  portion  of  this  investigation',  the 
instantaneous  velocity  components  in  the  wakes  of  the  65  AF 
and  120  AF  propellers  at  one  pitch  setting  each  were  mea¬ 
sured  by  tneans  of  hot-wire  anemometry.  Tne  Canadair,  Ltd., 
static  thrust  facility  was  used,  and  the  performance  wqs 
measured  for  these  propellers  over  a  range  of  pitch  settings. 
The  full  details  of  the  acquisition  and  reduction  of  the 
data  are  contained  in  the  test  report  prepared  by  D.  C. 

Gilmore  and  I.  S.  Gartshore  of  McGill  University,  Ref.  14. 

However,  subsequent  to  the  preparation  of  Ref.  14,  the 
basic  velocity- voltage  calibration  procedure  was  reexamined. 
This  led  to  a  redetermination  of  the  calibration  constants 
in  terms  of  true  voltage  instead  of  the  voltage  read  on  the 
hot-wire  control  unit.  The  velocity  components,  as  reduced 
on  the  basis  of  the  revised  calibration,  showed  appreciable 
changes,  which  will  be  seen  in  the  presentation  of  results 
that  follows.  These  changes  display  clearly  the  importance 
of  accurate  velocity- voltage  calibrations,  which  is  emphasized 
in  Ref.  14. 

The  instantaneous  axial  and  radial  components  of  the 
velocity  at  various  field  points  behind  the  65  AF  propeller, 
operating  at  the  8.2°  pitch  setting  and  a  tip  speed  of  800 
feet  per  second,  are  presented  in  Figure  19.  The  corre¬ 
sponding  results  for  the  120  AF  propeller  at  the  10.0°  set¬ 
ting  and  a  700- foot-per-second  tip  speed  are  given  in  Figure 
20.  These  components,  nondimensionalized  with  propeller  tip 
speed,  are  plotted  as  a  function  of  the  azimuthal  coordinate 
0  in  propeller- f ixed  coordinates,  rather  than  as  a  function 
of  time  in  the  space- fixed  coordinates  in  which  they  were 
actually  measured.  Note  that  for  convenience  the  negative 
of  the  radial  component  has  been  plotted.  The  tangential 
component  is  not  presented  because  it  is  so  much  smaller  in 
magnitude  than  the  other  two  components  that  it  cannot  be 
resolved  accurately;  see  Ref.  14.  The  data  that  are  pre¬ 
sented  were  reduced  on  the  basis  of  the  revised  calibration 
constants,  with  the  original  results  of  Ref.  14  included  for 
comparison.  The  algebraic  sign  of  the  differences  between 
results  based  on  the  original  calibrations  and  those  based 
on  the  revised  calibrations  depends  upon  the  particular  hot¬ 
wire  probe  used.  The  significance  of  the  repeated  data 
points  is  discussed  on  pp.  29  and  30  of  Ref.  14. 
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Trailing  Vortex  Sheet  Location 
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Azimuthal  Coordinate 
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Azimuthal  Coordinate 


x/Rp=0.1,  r/Rp=0 . 75  x/Rp=0.1,  r/Rp=1.0  x/Rp0.25,  r/R, 

Azimuthal  Coordinate,  0 


sjuauoduioo  X^tooisa  <atiU/A  FUB  <aHy/n  ‘TBixy 


Figure  20.  Contd. 


For  further  discussions  of  the  wake  data,  see  the 
sections  in  Chapters  4  and  6  concerning  the  65  AF  propeller 
at  the  8.2°  pitch  setting  and  the  120  AF  propeller  at  the 
10.0°  pitch  setting,  respectively. 
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